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BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 

During  the  last  six  month  period  of  the  Grant,  work  was  performed 
on  the  spectroscopy  of  water  in  both  liquid  and  vapor  state.  These 
investigations  were  performed  with  an  optically  pumped  laser  as  the 
submillimeter  radiation  source.  The  details  of  this  work  have  been 
presented  in  the  following  Conference  papers  and  journal  submissions: 


0.  A.  Simpson,  B.  L.  Bean  and  S.  Perkowitz,  "Far  Infrared  Optical 
Constants  of  Liquid  Water  Measured  With  An  Optically  Pumped  Laser," 
accepted  for  publication  in  the  Journal  of  the  Optical  Society  of 
America  -  December,  1979  issue. 

0.  A.  Simpson,  R.  A.  Bohlander,  J.  J.  Gallagher  and  S.  Perkowitz, 
"Measurements  of  Far  Infrared  Water  Vapor  Between  Lines  With  An 
Optically  Pumped  Laser,"  presented  at  The  Conference  on  Microwave 
Spectroscopy  and  Coherent  Radiation,  Duke  University,  June,  1979 
and  submitted  to  the  Journal  of  Physical  Chemistry. 

R.  A.  Bohlander  et  al ,  "Excess  Absorption  by  Water  Vapor  and 
Comparison  With  Theoretical  Dimer  Absorption,"  Proceedings  of  the 
Workshop  on  Water  Vapor,  Vail,  Colorado,  September  11-13,  1979. 

0.  A.  Simpson,  R.  A.  Bohlander,  J.  J.  Gallagher  and  S.  Perkowitz, 
"Submillimeter  Spectroscopy  of  Water  Vapor,”  Conference  on 
Microwave  Spectroscopy  and  Coherent  Radiation,  Duke  University, 
June,  1979. 

0.  A.  Simpson,  B.  L.  Bean  and  S.  Perkowitz,  "Submillimeter  Laser 
Spectroscopy  of  Liquid  Water,”  The  Conference  on  Microwave 
Spectroscopy  and  Coherent  Radiation,  Duke  University,  June,  1979. 

J.  J.  Gallagher,  "Applications  of  Millimeter  and  Submillimeter 
Spectroscopy  to  Atmospheric  Propagation,”  Conference  on  Microwave 
Spectroscopy  and  Coherent  Radiation,  Duke  University,  June,  1979 
(Invited  Paper). 

0.  A.  Simpson,  S.  Perkowitz,  R.  A.  Bohlander  and  J.  J.  Gallagher, 
"Far-Infrared  Laser  Spectroscopy  of  Water  Vapor  and  Liquid  Water," 
Digest  of  the  Fourth  International  Conference  on  Infrared  and 
Millimeter  Waves,  Miami,  Florida,  December  10-15,  1979. 
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These  measurements  with  the  optically  pumped  laser  are  the  initial 
measurements  of  this  nature  and  will  be  extended  to  obtain  detailed  data 
on  atmospheric  absorption.  The  work  is  part  of  the  requirements  for  the 
doctorate  degree  of  0.  A.  Simpson. 

A  quasi-optical  dlplexer  and  bi-conieal  Schottky  barrier  mixer  have 
been  constructed  and  are  undergoing  tests  with  an  optically  pumped  laser 
as  local  oscillator.  This  work  will  be  reported  when  final  results  are 
obtained . 

Two  manuscripts  currently  in  preparation  have  been  partially 
supported  by  ARO.  Copies  of  these  will  be  forwarded  to  the  ARO  Physics 
Division  upon  completion.  They  are  a  paper  on  Submillimeter  Wave 
Techniques,  for  the  Proceedings  of  the  IEEE,  and  one  on  Quasi-Optics, 
which  will  be  a  chapter  in  the  Academic  Press  series  on  Millimeter 
through  IR  Physics. 

During  the  course  of  this  program,  many  topics  on  submillimeter 
spectroscopy  and  techniques  have  been  investigated  and  reported  in 
previous  semi-annual  reports.  Many  of  these  topics  have  potential 
applications  in  the  area  of  NMMW  military  systems.  A  paper,  presented 
at  the  1979  IRIS  Conference  and  at  the  1979  SP1E  Conference  as  invited 
papers,  shows  the  importance  of  spectroscopy  related  to  atmospheric 
effects  and  research  on  technology  in  this  spectral  region.  This  paper 
is  a  result  of  partial  support  by  HDL  under  ARO  Contract  No.  DAAG29-77- 
C-0026.  A  technical  report  presenting  an  overview  of  the  research 
performed  under  this  grant  is  in  preparation  and  will  be  submitted  to 
the  ARO  Physics  Division. 


A 


ACCEPTED  FOR  PUBLICATION  IN  THE  JOURNAL 
OF  THE  OPTICAL  SOCIETY  OF  AMERICA 
FINAL  VERSION 

FAR  INFRARED  OPTICAL  CONSTANTS  OF  LIQUID  WATER 
MEASURED  WITH  AN  OPTICALLY  PUMPED  LASER* 


0.  A.  Simpson,  3.  L.  Bean,+  and  S.  Perkowitz 
Physics  Department 
Emory  University 
Atlanta,  Georgia  30322 


A  tunable  far  infrared  optically  pumped  laser  has  been  used  to  measure 
the  reflection  R  and  transmission  of  liquid  water  at  25°C.  The  laser 
covered  the  range  8.22  to  175-7  cm"1  with  an  average  spacing  between 
lines  of  8  cm-1.  The  optical  data  were  used  to  find  the  complex  index 
of  refraction  n-ik  with  typical  errors  of  2^.  A  small  peak  in  n  was 
observed  near  55  cm-1.  This  feature  correlates  with  published  Raman 
data.  The  values  of  R,  n,  k  and  the  absorption  coefficient  a  are  given 
in  graphical  and  tabular  form. 


*Work  partially  supported  by  the  Army  Research  Office  under  Grant  No. 
DAAG29-76-G-0280. 

tPresent  address:  Science  Applications  Inc.,  Drawer  E,  White  Sands  Missile 


Range,  New  Mexico  88002. 


T.  INTRODUCTION 


Liquid  water  absorbs  so  heavily  in  the  far  infrared  (FIR)  range  that 
it  :s  difficult  to  accurately  measure  its  transmission.  Many  workers 
have  determined  the  complex  index  of  refraction  n-ik  with  varying 
degrees  of  precision.^  Afsar  and  Hasted  ( AH )  have  avoided  the  absorp¬ 
tion  problem  by  using  reflection  dispersive  Fourier  transform  spectroscopy. 
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They  first  measured  n-ik  at  a  water  temperature  of  19°C  and  later  accurately 

found  n-ik  at  L,  30,  and  57°C.'5  An  alternate  method  is  to  directly  measure 

the  transmission  through  a  siiable  water  pathlengt’n  with  a  high  power 

optically  pumped  FIR  laser.  Since  free  surface  reflection  measurements 

are  also  feasible  with  a  laser,  FIR  laser  spectroscopy  gives  n-ik  in  an 

elegantly  simple  manner.  In  this  paper  we  present  laser  results  between 

8.22  and  175-7  cm-1  at  a  water  temperature  of  25°C.  Our  data  are  superior 

k 

to  the  previously  reported  results  at  this  temperature,  provide  better 
low  frequency  coverage  than  was  available  to  AH,  and  give  evidence  for  a 
temperature  dependent  low  frequency  feature  near  55  cm-1.  The  general  agreement 
between  the  results  of  AH  and  our  data  is  so  good  as  to  confirm  the  validity 
of  both  experimental  techniques. 


IT.  EXPERIMENTAL  METHODS  AND  RESULTS 


A  block  diagram  of  the  experimental  apparatus  is  shown  in  Fig.  1. 

A  25  W  CO^  laser  drives  a  FIR  laser  of  the  waveguide  type.  Details  of 
this  laser  spectrometer  (and  early  water  measurements)  have  been  previously 
reported,  but  some  improvements  have  been  made.  A  new  C0^  laser  was 
constructed  with  a  PZT  mounted  ZnSe  mirror  to  allow  tuning  of  the  line 


■*V<:  • 


■m  hh  ■■ 
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center  frequency  to  match  the  absorption  transition  frequency.  Sampling 

of  the  FIR  oeam  (signal  I  ) b  •  the  pyroelectric  detector  D  allowed  the 

s  1 

C0o  laser  to  be  frequency  stabilized  so  as  to  keep  the  FIR  power  constant. 

This  helped  to  eliminate  amplitude  fluctuations  in  the  FIR  laser  output 
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caused  by  instabilities  in  the  CO^  laser.  An  opto-acousti c  detector  on 
the  FIR  cavity  aided  in  tuning  difficult  absorption  transitions. 

The  main  signal  entered  the  Fabry-Perot  interferometer  (FP)  which 
was  used  both  to  measure  the  wavelength  and  as  a  filter  to  eliminate 
unwanted  wavelengths  that  may  have  lased  simultaneously  with  the  desired 
signal.  The  FIR  beam  was  then  divided  by  BS^.  The  reference  signal  1^, 
after  detection  and  amplification  by  a  Golay  cell  and  lock-in  amplifier, 
entered  the  denominator  channel  of  the  ratiometer  RA.  The  signal  I 
was,  as  will  be  described,  reflected  from  or  transmitted  through  the 
water  sample  and  then  detected  by  another  Golay  cell.  This  signal, 
after  lock-in  amplification,  became  the  numerator  signal  in  the  ratiometer. 
This  ratioing  process  was  important  in  eliminating  additional  source 
fluctuations  in  the  FIR  power  level.  Values  of  1^/1^  from  the  ratiometer 
were  the  basic  experimental  signal. 

Although  there  are  hundreds  of  optically  pumped  FIR  lines  arising 
from  dozens  of  pumped  media,  it  is  simplest  to  work  with  only  a  few  media 
when  possible.  We  obtained  23  FIR  lines  between  8.22  and  175*7  cm-1  using 

g 

only  four  pumped  gases  as  has  been  fully  described  elsewhere.  This 
frequency  range  and  coverage  is  comparable  to  that  available  to  AH 
in  their  dispersive  Fourier  work. 

Free  surface  measurements  of  the  reflection  coefficient  R  were  made 
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mirror  deflected  the  KIR  beam  to  near- 


ontal  water  surface  and  a  second  mirror 


directed  the  reflected  beam  into  the  detector 


was  snail  enough  that  the  assumption  of  normal  incidence  introduced  a 


R  radiation  was  unnolarized  which 


ible  error.  In  additi 


further  decreased  this  error.  The  water  was  placed  in  a  container 


P,  so  that  the  back  surface  reflection  was  totally 


negligible.  The  water  container  could  be  replaced  when  desired  with  a 


reflectivity.  The 


standard  reflector 


coincidence  of  the  water  and  mirror  surfaces  on  interchange  was  constantly 


monitored  by  the  reflection  of  a  fixed  HeNe  laser  beam 


For  each  FIR  frequency  several  values  of  I./IR  were  determined  for 


ce  and  reflection  from  the  standard  mirror 


reflection  from  the  water  sur 


The  ratio  of  the  two  values  of  I . / I  then  yielded  R  with  typical  random 


’able  I 


under  1 


The  measure5,  values  are 


An  adjustable  pathlength  cell  with  crystalline  quartz  windows  wa: 


The  cell  was  mounted  on  a  t ran s la- 


used  in  the  transmission  men 


tage  arrangement  with  a  high  accuracy  (0.002  mm/division) 


dial  gauge,  attached  directly  to  the  cell 


micrometer  head  screw  drive 


gave  a  second  determination  of  pathlength  changer,  for  comparison 


were  measured  at  increments 


with  the  micrometer  drive 


r  an  average  change  in  pathlength  of  60pm 


of  to  an  accuracy  of  1pm  o 


These  data  made  it  possible  to  determine  the  Lambert  absorption  coefficient 


a  from  the  relation 


where  x  is  the  puthlength.  An  error  weighted  least  squares  fit  was  applied 
to  the  data  at  each  laser  frequency  f  to  determine  a(f)  and  the  associated 
error,  typically  2 i*.  The  resulting  values  of  a  are  shown  in  Fig.  2. 

Values  of  k  =  a/Unf  were  calculated  from  the  measured  absorptions  and 
together  with  the  appropriate  R,  were  used  to  determine  n  from  the  normal 
incidence  expression: 


t<g|>  -  <K!*1)],/! 


(2) 


The  results  are  plotted  in  Fig.  3>  Numerical  values  of  all  the  measured 
and  derived  optical  quantities  a.te  given  in  Table  I. 


III.  DISCUSSION 

Our  data  at  25°C  are  in  excellent  agreement  with  those  of  AH'  at  30°C. 
Other  published  data  near  25°C  disagree  strongly  with  these  combined  results. 
The  19°C  values  AH  found  for  n  are  veil  below  the  present  data  between  50 
and  150  cm-1  and  even  lie  below  their  own  later  k°C  data  between  TO  and 
150  cm"  .  The  Downing  and  Williams^  values  for  n  at  27°C,  on  the  other 
hand,  are  significantly  greater  than  the  present  results  above  100  cm-1. 

These  conflicting  observations  could  be  reconciled  only  if  the  FIR  properties 
of  water  show  radical  and  non-monotonic  temperature  dependences  between  1*  and 
30°C.  This  possibility  remains  to  be  explored,  but  the  good  agreement  between 
the  present  data  at  25 °C  and  the  AH  data  at  30°C  suggests  that  the  results  at 
19  and  2T°C  are  experimental  anomalies. 

An  interesting  feature  of  the  laser  results  is  the  small  maximum  in  n 
near  55  cm-1.  Raman  measurements^*  ‘ ^  have  shown  the  presence  of  a  small 
peak  near  60  cm-1  at  25°C  but  earlier  FIR  data  do  not  exhibit  this  feature 
because  of  the  relatively  poor  signal-to-noise  ratio.  The  peak  magnitude  is 
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clearly  greater  than  the  raniom  error  in  our  measurements ,  however. 

IV.  C0?ICLU3I0NS 

Our  laser  measurements  s*  n  and  k  at  ?C<°C  are  in  excellent  agreement  with 
re:'  ■  Fourie  lata  of  A!  r  fttu!  Hss'.t'l.  ;  j  ■ 

error-.;  in  the  measurements  are  comparable.  Assuring  that  there  is  no  strong 
ter.; -erat  .re  dependence  in  the  genera1,  optical  properties  of  water  in  this 
tempo  rat sro  range,  th.e  agreement  attests  to  the  validity  of  both  experimental 
techniques . 

Our  data  exhibit  a  low  frequency  peak  near  55  cm- 1  which  appears  to 
correlate  with  similar  structure  seen  in  Raman  data.  The  Raman  peak  is 
believed  to  arise  from  hydro -on  bond  stretching  and  bending  notions.^*"  We 
have  begun  to  examine  the  peak  behavior  by  making  multi-oscillator,  relaxation 
model,  and  Kraners-Kronig  analyses  of  our  highly  accurate  reflectivity  data. 
Preliminary  results  suggest  that  relaxation  rather  than  resonant  effects 
dominate  in  this  frequency  region.  This  analysis  is  continuing  and  we  plan  also 
to  study  the  position  and  amplitude  of  th.e  peak  vs.  temperature.  The  laser 
spectrometer  will  be  used  for  further  accurate  measurements  as  a  function  of 
temperature, with  potential  for  shedding  light  on  basic  molecular  processes 
in  water. 
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'fjy  vi'i 


f(cm  *) 

X(vim) 

R 

,  -lx 
a  (cm  ) 

n 

k 

8.2 

1217 

0.297 

121 

2.857 

1.174 

11.2 

890 

0.249 

150 

2.462 

1.062 

13.1 

764.1 

0.228 

148 

2.437 

0.902 

14.3 

699.5 

0.228 

171 

2. 383 

0.950 

17.5 

570.5 

0.201 

173 

2.  306 

0.784 

21.2 

471 

0.191 

199 

2.252 

0.746 

24.0 

417 

0.185 

215 

2.231 

0.714 

24.6 

406.4 

0. 184 

206 

2.267 

0 . 665 

26.8 

372.7 

0.179 

235 

2.197 

0.696 

30.3 

330 

0.164 

228 

2.155 

0.598 

34.2 

292.6 

0.161 

259 

2.126 

0.604 

35.5 

281.5 

0.155 

250 

2.112 

0.560 

39.4 

254 

0.147 

221 

2.125 

0.447 

43.6 

229.1 

0.144 

289 

2.050 

0.527 

52.1 

191.7 

0.140 

285 

2.079 

0.435 

54.5 

183.6 

0.147 

355 

2.082 

0.519 

61.3 

163 

0.140 

339 

2.076 

0.440 

69.1 

144.8 

0.136 

405 

2.025 

0.466 

84.2 

118.8 

0.134 

457 

2.034 

0.432 

98.4 

101.6 

0.131 

605 

1.976 

0.489 

103.6 

96.5 

0.136 

696 

1.985 

0.534 

141.6 

70.6 

0.122 

986 

1.845 

0.554 

175.7 

56.9 

0.095 

1147 

• 

J ■ 
N) 

0.519 

Tabic  I.  Measured  values  of  reflection  R  and  absorption  coefficient  a,  and 

derived  values  of  the  Index  of  refraction  n  Ik  for  liquid  water  at  ?s*C. 
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FMURF  CAPTIONS 


Fi*urp  1  '  Pl0ck  diagram  of  the  FIR  laser  spectrometer t  m,  mirrors;  0, 
^ro*i„  ;  Una;  SA,  sreet.ru-  analyser;  PST,  piezoelectric  transducer; 

L*  0,‘°pp<*r»  Jo  tec  tor*;  RP,  beamsplitters;  FP,  Fabry -Perot  interferom- 
oter,  s,  water  sample;  1,1,  look-in  amplifier;  HA,  digital  rat  low?  ter;  FP , 

frequency  stabiliser;  PC,  stripchart  recorder;  OAD,  opto-aeoustic  detec- 
to  r . 


hpU,v  T,a"r>r  t'eaults  for  the  absorption  coefficient  of  liquid  water 
at  ?5°C.  The  symbols  for  the  laser  experimental  points  are  slightly 
larger  than  the  typical  ran  do-  errors. 

Hgure  3.  Laser  results  for  the  real  part  of  the  index  of  refraction  of 
liquid  water  at  Th°C.  The  symbols  for  the  laser  experimental  points 
are  slightly  larger  than  the  typical  random  errors. 
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ABSORPTION  BETWEEN  LINKS  WITH  AN  OPTICALLY  PUMPED  LASER* 
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Georgia  Institute  of  Technology 
Atlanta,  Georgia  30332 

S.  Perkowitz 
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A  tunable  far-infrared  optically  pumped  laser  has  been  used  to  measure 
the  absorption  by  equilibrium  samples  of  water  vapor  at  25°C  in  the  line  gaps 
between  8.17  and  103.0  cm  The  experimental  results  are  compared  with 
calculations  of  predicted  absorption  based  on  the  line  shape  formula  given  by 
Gross.  Excess  absorption  lias  been  found  which  is  in  reasonable  agreement  with 
other  published  data.  New  data  are  presented  for  the  region  60  to  105  cm  ^ . 


*Work  supported  by  the  Army  Research  Office  under  Grant  No.  DAAG29-76-G-0280 
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I.  INTRODUCTION 

In  the  far  infrared  (KIR)  or  submillimeter  spectral  region  water 

vapor  is  the  dominant  absorber  in  the  atmosphere.  Previous  studies  of 

water  vapor,  both  in  the  laboratory  and  atmosphere  have  shown  absorption 

between  lines  which  is  in  excess  of  that  predicted  by  existing  line  shape 
12  3 

formulae.  '  ’  This  excess  absorption  has  been  attributed  to  errors  in 

1  2-7 

theoretical  line  shapes,  to  water  dimers,  and/or  molecular  aggregates 

g 

or  clusters.  Considering  the  recent  growth  of  atmospheric  applications 

in  the  FIR-submillimeter  region,  the  lack  of  a  firm  understanding  of  how 

radiation  and  water  molecules  interact  only  underscores  the  need  for  more 

hard  data  from  controlled  experiments. 

The  present  study  is  an  initial  effort  to  measure  the  absorption  by 

equilibrium  samples  of  water  vapor  of  radiation  from  a  tunable  FIR  optically 
9 

pumped  laser.  The  use  of  a  pumped  FIR  laser  is  a  unique  feature  of  our 
measurements.  The  high  power  of  this  device  gives  a  signal-to-noise  ratio 
which  is  superior  to  any  other  available  broadband  FIR  source. 

Data  are  presented  for  various  water  vapor  pressures  at  23°C  for  seven 
laser  frequencies  between  8.17  and  103.6  cm  The  frequencies  chosen  lie 
in  the  gaps  between  strong  monomer  lines  as  shown  in  Fig.  1  and  were  taken 
from  the  National  Physical  Laboratory  table. ^  Only  those  frequencies  well 
removed  from  line  shoulders  were  chosen  in  order  to  minimize  the  error  in 
predicted  absorption  due  to  uncertainty  in  the  laser  frequency.  Although 
the  family  of  discrete  FIR  laser  lines  does  not  provide  continuous  spectral 
coverage,  there  was  no  difficulty  in  finding  lines  which  met  this  criteria. 
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Our  measured  absorptions  are  In  excess  of  the  theoretical  values  at 
each  of  the  laser  frequencies  and  at  all  pressures.  The  derived  excess 
absorpt ion  agrees  reasonably  well  with  other  published  data  obtained  with 
Fourier  spectrometers  and  HCN  lasers  andorovldes  new  information  in  the 
spectral  region  between  00  and  105  cm  *. 


II.  THEORETICAL  AND  EXPERIMENTAL  METHODS 


The  predicted  absorption  at  each  laser  frequency  was  calculated  using 

11  12 

the  line  shape  formula  given  by  Gross  ’  : 


,  ^  nS  4v  6 

a(v)  =  —  — ^ - — 

(v  -  v  )  +  4v  5 

u  o 


for  the  contribution  of  a  line  centered  at  wavenumber  vq  to  the  absorption 
coefficient  a  at  wavenumber  v.  Here  <5  is  the  half-width  at  half  maximum,  S 
is  the  line  intensity,  and  n  is  the  absorbing  medium  number  density.  There 
is  very  little  difference  in  choosing  the  Gross  formula  over  other  commonly 
used  line  shape  formulae  (Van  Vleck-Weisskopf ,  Lorentz,  etc.).  These  are 
all  based  on  instantaneous  impact  theory  and  give  essentially  idential  result 
in  the  spectral  region  studied.  The  formula  given  by  Van  Vleck  and  Weisskopf 
does  predict  less  absorption  below  15  cm  1,  but  in  this  region  excess  absorp¬ 


tion  is  large  enough  that  the  difference  is  not  significant. 

The  water  line  frequencies,  intensities,  and  initial  state  energies  for 
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the  346  lines  used  were  taken  from  the  AFGL  table  with  self-broadening 
line  width  coefficients  taken  from  calculations  by  Benedict  and  Kaplan^. 

The  modification  of  line  parameters  for  temperature  and  pressure  was  done 
as  described  by  McClatchev  with  the  exception  of  the  temperature  dependent 
exponent  of  the  line  widths,  which  were  determined  separately  for  each 
line  from  values  of  widths  tabulated  at  two  temperatures  by  Benedict  and 
Kaplan. 

The  measurements  were  made  with  an  optically  pumped  waveguide  type 
FIR  laser  spectrometer,  which  is  shown  in  block  diagram  form  in  Fig.  2. 

This  apparatus  has  been  e  xtensively  described  elsewhere***’  ’  ***  and  only 
relevant  details  will  be  given  here.  The  laser  provides  a  large  selection 
of  individual  FIR  lines  with  a  typical  power  of  milliwatts  CW  in  each  line. 

The  frequencies  have  been  determined  with  great  accuracy  by  various  workers, 
but  as  a  precaution  against  possible  error  in  the  published  laser  frequencies, 
a  Fabry-Perot  interferometer  was  used  to  measure  and  confirm  each  of  the 
seven  wavelengths. 

The  water  vapor  was  contained  in  a  straight  glass  pipe  absorption  cell 
of  length  3.44  m  and  diameter  10  cm.  The  ends  of  the  cell  had  crystalline 
quartz  windows  offset  15°  to  the  laser  beam  path.  Initial  measurements  with 
the  cell  had  resulted  in  much  larger  absorption  than  expected,  which  was 
determined  to  be  due  to  interference  effects.  The  offset  of  the  windows  was 
helpful  in  eliminating  these  interference  effects,  as  were  several  co-axial 
absorbing  bellows  placed  at  several  points  along  the  cell.  Distilled  water 
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was  evaporated  and  introduced  into  the  cell  through  a  perforated  tube 
running  the  cell  length.  When  sufficient  vapor  had  entered  the  cell,  the 
system  was  closed  and  allowed  to  reach  equilibrium  at  room  temperature, 

25  +  1°C.  The  static  vapor  pressure  was  measured  with  a  Baratron  MKS 
capacitance  manometer,  which  had  an  accuracy  of  +2%. 

To  measure  transmission  the  output  from  a  Golay  detector  preceding 
the  absorption  cell  was  electronically  divided  into  the  output  of  another 
Golay  detector  following  the  cell.  For  each  FIR  laser  frequency,  several 
values  of  this  ratio  were  determined  for  transmission  through  the  cell 
evacuated  and  for  various  water  vapor  pressures.  The  resulting  ratio  of 
the  two  separately  determined  values  then  yielded  the  measured  water  vapor 
power  transmission  coefficient.  This  method  has  the  important  advantage 
of  eliminating  any  fluctuations  in  the  FIR  power  arising  from  laser 
instabilities.  Measured  values  of  the  water  vapor  absorption  at  A3. 649  cm 
are  shown  in  Fig.  3  as  an  example,  together  with  the  predicted  absorption 
based  on  the  Gross  line  shape. 


The  excess  absorption  coefficient  was  calculated  according  to  the  equation 

a  -  -en(TM/T„)  /  L 

excess  M  G 

where  L  is  the  pathlength,  T  and  T^  are  the  measured  and  theoretical 
transmission  respectively.  As  shown  in  Fig.  4  for  the  frequency  43.649  cm  ^ , 
a  straight  line  showing  excess  absorption  depending  on  density  squared  fits 


the  data  within  experimental  undertainty.  Such  a  relationship,  determined  by 


an  error  weighted  least  squares  fit  was  found  for  each  of  the  other  laser 
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frequencies.  This  result  was  expected,  since  in  both  speed  of  sound  and 
20 

calorimetric  data  the  third  and  higher  order  virial  terms  have  been 
reported  to  be  negligible.  The  observed  excess  absorption,  normalized  to 
density  squared,  is  shown  in  Fig.  5  for  each  of  the  laser  frequencies. 

III.  DISCUSSION 

In  a  study  of  molecular  absorption  one  has  to  keep  in  mind  the 

limitations  of  present  line  shape  formulae  that  are  based  on  the  assumption 

of  instantaneous  collisions.  Excluded  from  these,  of  course,  are  a  wide 

range  of  possible  molecular  interactions.  This  particular  weakness  in  theory 

introduces  difficulties  in  interpreting  new  information;  however,  such 

formulae  are  widely  used.  Hence,  calculations  in  the  present  study  can  be 

readily  reproduced  by  other  researchers. 

Our  measurements  of  the  excess  absorption  in  the  FIR  spectral  region 

by  water  vapor  are  in  general  agreement  with  other  reported  work,  but  we 

have  observed  a  somewhat  larger  excess  absorption  at  50  cm  It  should 

2 

be  noted  that  similar  studies  on  water  vapor  with  a  HCN  maser  also 
resulted  in  slightly  greater  excess  absorption  that  that  from  broadband 
FIR  source  measurements.  In  addition,  there  is  new  evidence  of  an  unexpected 
decrease  in  the  excess  absorption  at  higher  frequencies.  These  results  are 
preliminary  and  more  laser  lines  will  be  used  to  increase  resolution  in 
this  frequency  region. 
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FIGURE  CAPTIONS 


Figure  la,  b.  Theoretical  water  vapor  absorption  based  on  the  Gross 
line  shape  formula  with  laser  frequencies  in  the  gaps 
shown  by  arrows. 

Figure  2.  Block  diagram  of  the  FIR  laser  spectrometer:  m,  mirrors; 

G,  grating;  L,  lens;  SA,  spectrum  analyzer;  PZT  piezo¬ 
electric  transducer;  C,  chopper;  D,  detectors;  BS,  beam¬ 
splitters;  FP,  Fabry-Perot  interferometer;  S,  water  vapor 
sample;  LI,  lock-in  amplifier;  RA,  digital  ratiometer; 

FS,  frequency  stabilizer,  SC,  stripchart  recorder;  OAD, 
opto-acoustic  detector. 

Figure  3.  The  measured  water  vapor  absorption  at  43.649  cm  ^ . 

Theoretical  absorption  shown  by  X's  for  comparison. 

Figure  4.  Density  dependence  of  excess  absorption  at  43.649  cm 

Straight  lines  fitted  by  weighted  least  squares. 

Figure  5.  Observed  excess  absorption  coefficients  divided  by  den¬ 

sity  squared  for  each  laser  frequency. 
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The  predicted  absorption  was  calculated  at 
each  laser  frequency  vising  the  line  shape 
formula  given  by  Gross  [4,5]: 


Introduction 

Wo  have  used  a  far  Infrared  (KIR)  optically 
pumped  laser  to  measure  both  the  reflection  and 
transmission  of  liquid  water  and  the  absorption 
by  water  vapor  of  radiation  In  the  spectral  re¬ 
gion  betwen  8.22  and  175.7  cm”1.  The  pumped 
laser  is  a  unique  feature  of  our  measurements. 

The  high  power  of  this  device  gives  a  signal- to* 
noise  ratio  which  Is  superior  to  that  available 
from  broadband  FIR  sources. 

The  optical  data  for  liquid  water  were  used 
to  determine  the  complex  refractive  index 
n-lk  [ 1  ] .  Data  for  measured  absorption  by  water 
vapor  were  compared  with  calculations  of  pre¬ 
dicted  absorption  based  on  the  line  shape 
formula  given  bv  Gross.  Our  measured  absorptions 
are  In  excess  of  the  theoretical  values  at  each 
of  the  laser  frequencies  and  at  all  pressures. 


for  the  contribution  of  a  line  centered  at 
frequency  vc  to  the  absorption  coefficient  K  at 
frequency  v.  Here  6  is  the  half-width  at  half 
maximum,  s  Is  the  line  Intensity,  and  n  is  the 
absorbing  medium  number  density.  The  water  line 
frequencies,  intensities,  and  initial  state 
energies  were  taken  from  the  AFGL  [b]  table  with 
self-broadening  line  width  coefficients  taken 
from  calculations  by  Benedict  and  Kaplan  [7]. 

The  modi f icat ion  of  line  parameter-  for  tempera¬ 
ture  and  pressure  was  done  as  described  in  the 
AFGL  table  with  the  exception  of  the  temperature 
dependent  exponent  of  the  line  widths,  which 
w**-e  determined  separately  for  each  line  from 
values  of  width  tabulated  at  two  temperatures  bv 
Benedict  and  Kaplan. 

The  excess  absorption  coefficient  was 
calculated  according  to  the  equation 


Experimental  Methods  and  Results 

The  FIR  laser  was  of  the  waveguide  tvpe  and 
has  been  described  elsewhere  [1-3].  A  ratioing 
process  using  two  detectors  was  helpful  in 
eliminating  anv  fluctuations  In  the  FIR  power 
arising  from  laser  Instabilities. 

Free  surface  measurements  of  the  reflec¬ 
tion  coefficient  R  were  made  on  distilled  water 
at  25  ♦  1°C  for  near  normal  incidence  with 
typical  random  errors  under  IX.  An  Adjustable 
path  length  cell  was  used  to  determine  values  of 
the  liquid  WAter  transmission  at  increments  of 
5  urn  to  an  accuracy  of  l  um  over  an  average 
change  in  pathlength  of  60  urn.  The  Lambert 
absorption  coefficient  o  was  determined  at  each 
laser  frequency  f  by  least  squares  fittings  and 
used  to  calculate  values  of  k  -  o/4nf.  Each 
k  with  the  appropriate  R  was  used  to  determine 
n  from  the  normal  Incidence  relation.  The 
results  for  n  are  shown  In  Fig.  1. 

Measurements  of  the  absorption  hy  water 
vapor  were  performed  using  a  straight  glass 
pipe  absorption  cell  of  length  3.44  m  and 
diameter  10  cm.  Distilled  water  was  evaporated 
and  introduced  into  the  cell  through  a  perfo¬ 
rated  tube  running  the*  cell  length.  To  measure 
transmission  the  output  from  a  Go lay  detector 
preceeding  the  absorption  cell  was  electronically 
divided  into  the  output  of  another  Go lay  detec¬ 
tor  following  the  cell.  The  quotient  of  the  two 
ratios  for  transmission  through  the  evacuated 
cell  and  for  various  water  vapor  pressures  yield¬ 
ed  the  measured  water  vapor  power  transmission 
coefficient . 


excess  H  G 

where  L  is  the  pathlength  and  T„  and  T^  are  the 
measured  and  theoretical  transmission  respective 
ly.  Measured  values  of  the  water  vapor  absorp¬ 
tion  at  43.649  cm”1  are  shown  in  Fig.  2  as  an 
example,  together  with  the  predicted  absorption 
based  on  the  Gross  line  shape.  The  observed 
excess  absorption,  normalized  by  the  density 
squared  is  shown  in  Fig.  3  for  each  of  the  laser 
f requenc  tes . 


Conclusions 

The  FIR  laser  has  made  it  possible  to  obtain 
highly  precise  data  in  liquid  water.  The  liquid 
data  cover  a  broader  range  than  the  results  of 
Afsar  [8],  hut  agree  well  where  they  overlap. 
Those  two  sets  of  laser  data  are  the  definitive 
FIR  results  on  liquid  water  to  date. 

The  derived  excess  absorptions  of  water 
vapor  extend  the  frequency  range  of  other 
published  data  obtained  with  Fourier  spectrom¬ 
eters  ( 9 , 1 0 1  and  HCN  lasers  [ 10,11 ]  and  provide 
new  information  in  the  spectral  region  between 
60  and  105  cm”1. 

♦This  work  was  partially  supported  by  the  Army 
Research  Office. 
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Fig.  2.  The  measured  water  vapor  absorption 
at  43.649  cm-1.  Theoretical  absorption  shown 
by  X's  for  comparison. 


Fig.  3.  Observed  excess  absorption  coefficients 
divided  by  density  squared  for  each  laser 
frequency. 
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Fig.  1.  Laser  results  for  the  real  part  of  the 
Index  of  refraction  of  liquid  water  at  25°C.  The 
symbols  for  the  laser  experimental  points  .ire 
slightly  larger  than  the  typical  random  errors. 
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abstract 


Absorption  by  water  vapor  ha.  been  measured  l„  wavenumber  ranges  from 
8  105  cm  and  300  to  600  cm'1.  The  egress  over  prediction  with  the  Oroas 

line  shape  has  been  compared  with  theoretical  dimer  absorption.  The  tempera¬ 
ture  dependence  and  over-all  magnitude  of  excess  absorption  are  In  reasonable 
agreement  with  prediction  for  dimers,  but  the  shape  Is  different.  The  limita¬ 
tions  of  the  present  theory  of  dimer  absorption  are  discussed. 
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1.  Introduction 


It  Is  well-known  that  established  line  shape  formulas  based  on  an  Impact 
approximation  are  unable  to  predict  molecular  absorption  In  the  wings  of  ab¬ 
sorption  lines.  In  the  case  of  water,  observed  absorption  has  been  found  to 
exceed  prediction  in  the  gaps  between  lines  frcm  the  microwave  to  the  infrared 
regions  of  the  spectrum.  The  desire  for  a  better  theoretical  model  has  been 
spurred  by  interest  in  the  propagation  of  radiation  through  water  vapor  in  the 
atmosphere. 

To  get  an  Improved  theory  of  the  shape  of  water  vapor  absorption  lines, 
more  detailed  Information  about  molecular  interactions  during  collisions  will 
be  needed  than  is  available  at  the  present  time.  Since  two  water  molecules 
can  form  a  weak  hydrogen  bond,  dimers  have  also  been  suggested  (Viktorova  and 
Zhevakin  1967,  1971,  1975)  as  another  possible  source  of  excess  absorption. 
This  paper  reports  studies  of  excess  absorption  in  the  laboratory  in  the 
wavenumber  ranges  8  to  105  cm  *  and  300  to  600  cm  \  and  comparisons  are  made 
with  a  calculated  theoretical  spectrum  of  water  dimers.  This  is  in  the 
nature  of  a  progress  report  and  is  not  a  definitive  test  of  whether  dimers 
are  important  since  there  are  significant  limitations  in  our  understanding  of 
dimer  structure  at  normal  temperatures.  Many  of  the  details  of  this  study 
have  been  given  in  the  thesis  of  Bohlander  (1979)  and  further  publications 
will  follow. 

2.  Excess  Absorption — Definition 

One  cannot  measure  excess  absorption  per  se:  it  is  the  difference 
between  observed  absorption  and  prediction  based  on  a  line  shape  formula. 

Thus  it  is  important  to  define  carefully  how  the  prediction  is  made.  We 
have  used  the  parameters  for  the  frequency,  intensity,  and  widths  of  lines 
given  by  McClatchey  et  al.  (1973),  Benedict  and  Kaplan  (1963),  Flaud, 
Camy-Peyret  and  co-workers  (1976,  1977),  and  Toth  and  Margolis  (1975). 

Commonly  used  line  shape  formulas  are  the  Gross  or  kinetic  line  shape  formula 
(Gross  1955),  the  Van  Vleck-Weisskopf  formula,  and  the  modified  Lorentz 
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formula  (Van  Vleck  and  Weisskopf  1945),  and  from  these  there  Is  little  to 
choose  theoretically.  They  all  Invoke  the  key  assumption  that  collisions  are 
Instantaneous.  We  have  chosen  the  Gross  formula.  Numerical  results  with  the 
formulas  mentioned  are  negligibly  different  in  the  wavenumber  range  studied 
except  with  the  Van  Vleck-Weisskopf  formula  below  15  cm  It  predicts 
somewhat  less  absorption,  but  the  conclusions  of  this  paper  would  not  be 
changed  if  this  formula  had  been  used. 

Preliminary  predictions  of  the  absorption  spectrum  were  made  including 
all  the  pure  rotation  band  lines  in  the  calculation.  We  then  made  shortened 

tables  of  lines  from  which  predictions  in  50  cm  *  intervals  differed  negligi- 

* 

bly  from  a  full  calculation.  In  this  way  typically  400  and  no  more  than  700 
lines  were  needed  at  any  given  frequency  for  subsequent  calculations,  com¬ 
pared  with  2600  pure  rotation  lines  below  700  cm  ^  in  the  complete  tabulation 
we  used. 

Some  workers  have  used  much  more  abbreviated  lists  of  lines  in  their 
calculation  of  the  pure  rotation  line  absorption  (see  eg.  Burch  1968;  Burch 
et  al.  1974;  Frenkel  and  Wood  1964;  Gaut  and  Reifenstein  1971):  they  have 
selected  only  lines  near  the  frequency  of  interest.  Excess  absorption  in 
this  case  varies  slowly  with  frequency  and  can  be  incorporated  in  an  easy-to- 
use  empirical  formula  to  predict  absorption.  However,  in  physical  terms, 
these  selections  of  lines  are  rather  arbitrary  and  cause  confusion  when 
values  of  excess  absorption  are  compared,  since  allowance  must  be  made  for  the 
various  omissions  of  absorption  lines.  Moreover,  observations  have  shown 
that  some  of  the  empirical  formulations  have  oversimplified  the  dependence 
of  absorption  on  the  water  vapor  density  and  temperature.  We  have  determined 
these  dependences  from  observations  and  will  consider  them  in  relation  to  the 
dimer  model. 


_3 

*The  tolerance  was  0.2  dB/km  for  a  water  vapor  density  of  26  g  m 


3.  Theoretical  Dimer  Absorption 


Ue  have  made  theoretical  calculations  of  the  absorption  spectrum  of  dimers 
from  existing  information  about  their  structure.  A  number  of  theoretical 
chemists  have  made  molecular  orbital  calculations  in  order  to  find  the  hydrogen 
bond  energy  and  the  most  stable  structure,  shown  in  Figure  1.  Dill  et  al. 

(1975)  give  a  useful  review,  and  more  recent  work  is  described  by  Matsuoka 
et  al.  (1976).  These  predictions  were  found  to  be  consistent  with  high  res¬ 
olution  measurements  (Dyke  et  al.,  1977)  of  microwave  transitions  between 
rotation  levels  in  the  ground  vibrational  state  of  dimers  in  molecular  beams. 
Therefore,  the  ground  vibrational  state  appears  to  be  well-understood  in  terms 
of  the  predicted  lowest  energy  structure. 

Since  the  hydrogen  bond  is  weak,  the  water  dimer  has  low-lying  vibration 

levels,  and  only  four  percent  of  the  population  is  estimated  to  be  in  the 

ground  vibrational  state  at  normal  temperatures.  Thus  it  is  an  important 

question  whether  dimers  retain  approximately  the  lowest  energy  configuration 

* 

at  these  temperatures.  Using  an  empirical  model  of  the  potential  energy  be¬ 
tween  two  water  molecules,  Owlcki,  Shipman  and  Scheraga  (1975)  have  concluded 
that  barriers  to  internal  rotation  are  low.  In  this  situation  the  orientations 
of  the  two  water  molecules  in  a  dimer  would  not  be  well-confined  to  configura¬ 
tions  near  the  lowest  energy  one.  We  felt  it  was  probably  premature,  and  was 
in  any  case  a  difficult  problem,  to  try  to  calculate  a  theoretical  dimer  spec¬ 
trum  that  would  take  account  of  internal  rotation.  Rather  we  have  calculated 
a  theoretical  spectrum  for  the  more  tractable  case  in  which  the  dimer  is  assumed 
to  undergo  rigid  rotation  and  small-amplitude  harmonic  vibration.  Comparison 
with  the  spectrum  of  excess  absorption  in  water  vapor  has  been  used  as  a  test 
of  this  theory  and  as  an  Indication  whether  there  may  be  the  need  to  consider 
internal  rotation. 

The  dimer’s  expected  pure  rotation  band  in  the  rigid  rotor  approximation 
comes  at  a  low  frequency  because  the  dipole  moment  (2.6  D)  measured  by  Dyke 
et  al.  (1977)  was  found  to  be  nearly  parallel  to  the  line  between  the  heavy 
oxygen  atoms.  The  peak  of  the  predicted  band  is  at  16  cm  (480  GHz)  as  shown 
in  Figure  2.  One  might  expect  the  spectrum  to  be  almost  that  of  a  symmetric 
top  with  a  rotation  inertial  constant  of  0.2  cm  ^  and  with  only  small  pertur¬ 
bations  due  to  the  light  hydrogen  atoms.  Early  calculations  of  the  dimer 
spectrum  were  made  on  that  basis  (Viktorova  and  Zhevakln  1967,  1971,  1975; 

£ 

Electrons  and  nuclei  were  placed  in  fixed  positions  to  reproduce  the  dipole 
moment  of  water  and  the  structure  of  ice. 
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Braun  and  Leidecker  1974).  However,  molecular  beam  studies  (Dyke  et  al.  1977) 
have  shown  a  much  more  complicated  spectrum  below  1.6  cm~*(48  GHz).  This  is 
attributed  to  inversion  tunnelling,  but  few  assignments  of  the  observed  lines 
have  been  made.  Therefore,  we  have  not  attempted  to  calculate  rotational 
structure  in  detail  and  have  used  smoothed  band  contours. 

The  remainder  of  the  predicted  absorption  is  assigned  to  six  intermolecu- 
lar  vibration  modes  involving  the  hydrogen  bond.  We  have  calculated  the  normal 
intermolecular  vibrational  frequencies  for  the  dimer  from  a  recent  molecular 
orbital  study  of  intermolecular  potential  energy  (Matsuoka  et  al.,  1976).  *  ** 
Most  of  the  modes  involve  partial  rotation  of  the  monomers  within  the  dimer, 
and,  since  this  causes  large  dipole  moment  oscillations,  the  predicted  absorp¬ 
tion  intensity  is  large.  Intensities  were  calculated  with  the  assumption  that 

the  charge  distribution  remains  fixed  on  the  monomer  units  during  vibra- 
** 

tions. 

In  addition,  we  need  to  know  the  number  density  of  dimers  to  be  able  to 
calculate  their  absorption  in  dB/km.  Although  dimers  have  been  observed  in 
molecular  beams,  we  know  nothing  from  this  about  equilibrium  concentrations. 
These  we  have  estimated  in  the  usual  way  (see  e.g.  Bolander  et  al.  1969)  from 
experimental  values  of  the  second  virial  coefficient  (Goff  and  Gratch,  1946). 
Typical  calculated  dimer  concentrations  are  of  the  order  of  a  part  in  1000  of 
the  monomer  concentration.  Uncertainties  in  the  theory  on  which  these  calcula¬ 
tions  are  based  will  be  discussed  later.  Dimer  concentrations  are  propor¬ 
tional  to  the  square  of  the  water  vapor  density,  and  their  dependence  on 
temperature  is  determined  by  the  hydrogen  bond  energy.  The  best  theoretical 
values  of  the  latter  are  around  .15  to  .16  eV;  the  value  estimated  from  the 
temperature  dependence  of  the  second  virial  coefficient  is  0.12  eV  (Bohlander, 
1979). 

* 

There  is  fairly  good  agreement  between  the  vibration  force  constants  from 
this  study  and  those  from  other  recent  molecular  orbital  studies  (Curtiss  and 

Pople  1975;  Klstenmacher  et  al.  1974). 

** 

Values  of  the  frequencies  and  intensities  are  also  in  reasonable  agreement 

with  those  found  by  (Xficki  et  al.  (1975)  using  an  empirical  potential  energy 
model. 


>.  Observations  In  the  Laboratory 

The  spectral  ranges  that  are  most  accessible  to  experimental  study  are  In 
the  wings  of  the  monomer  pure  rotation  band  and  Include  the  predicted  dimer 
pure  rotation  band  and  the  predicted  bands  for  hydrogen  bond  bending. 

We  have  studied  water  vapor  without  a  foreign  broadening  gas  since  this 
would  give  the  widest  possible  spectral  ranges  In  which  to  search  for  dimer 
absorption.  Uncertainties  in  the  theory  of  collision  broadening  by  a  foreign 
gas  are  also  avoided. 

We  began  the  study  using  Fourier  transform  spectroscopy  and  wlde-band 
radiation  taken  from  Mercury  lamp  or  glow  bar  sources.  Helium-cooled  bolo¬ 
meters  and  Golay  cell  detectors  were  used.  For  the  spectral  ranges  12  to  50 
cm  *  and  300  to  600  cm  *,  White-type  absorption  cells  provided  path  lengths 
between  20  and  200  meters.  A  selection  of  observed  spectra  is  given  in  Fig¬ 
ure  3.  They  represent  one  minus  the  ratio  of  spectra  obtained  with  a  sample 
and  with  a  vacuum  in  the  absorption  cell.  The  dashed  curves  show  expected 
monomer  absorption.  One  can  see  that  in  the  gaps  between  lines  in  these 
wavenumber  ranges,  excess  absorption  is  a  considerable  fraction  of  the  total. 
Whether  the  excess  absorption  was  due  to  collisions  between  two  molecules  or 
due  to  dimers,  it  was  expected  that  it  would  be  proportional  to  the  density 
squared,  and  this  was  found  to  be  true  within  experimental  certainty. 

Due  to  the  presence  of  numerous  strong  monomer  lines,  one  can  get  only  a 
fragmentary  picture  of  the  excess  absorption's  spectral  shape,  as  shown  for 
the  near  millimeter  region  in  Figure  A.  Some  results  of  Burch  (1968)  and 
Burroughs  et  al.  (1969)  are  shown  to  be  in  good  agreement.  Measurements 
were  also  made  by  one  of  us  (DU)  with  an  HCN  maser  and  are  shown  in  the 
figure  at  29.7  cm  Some  of  these  measurements  were  made  with  a  large 
untuned  cavity  which  Llewellyn-Jones  will  discuss  elsewhere  in  these 
proceedings.  He  will  also  present  further  data  in  the  range  2  to  15  cm-1 
obtained  with  this  facility. 

At  Georgia  Tech  and  Emory  University,  two  of  us  (0.  A.  S.  and  S.  P.) 
have  embarked  on  new  measurements  with  an  optically  pumped  laser  described 
elsewhere  (Bean  and  Parkowitx  1976,  1977).  Preliminary  results  are  shown 
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with  solid  circles.  There  is  the  potential  for  more  precise  measurements  be¬ 
cause  of  the  power  and  stability  of  the  laser,  and,  as  shown,  this  is  being 
used  to  extend  the  frequency  range  covered.  The  absorption  cell  is  a  straight 
pipe  3.4  m  long  designed  to  be  used  with  a  single  pass  through  it.  However, 
there  has  been  difficulty  in  eliminating  multiple  passes  arising  from  stray 
reflections.  This  has  been  largely  corrected  with  baffles,  but  there  is  pos¬ 
sibly  a  residual  effect  which  may  explain  why  these  data  show  somewhat  larger 
excess  absorption  at  lower  wavenumbers  than  previous  data. 

The  integral  of  the  observed  excess  absorption  shown  in  Figure  4  is  about 
the  same  as  that  predicted  by  the  dimer  model,  but  the  shape  of  the  spectrum 
is  different.  This  is  also  true  for  the  higher  wavenumber  range  studied 
(Figure  5) .  Here  observed  excess  absorption  is  represented  by  a  smooth  curve 
that  fits  our  data  and  that  of  Burch  et  al.  (1974)  within  the  experimental 
uncertainty.  As  previously  noted  by  Roberts  et  al.  (1976),  this  seems  to 
have  the  form  of  an  exponential  decrease  with  frequency  in  the  range  300  to 
600  cm  Further  discussion  of  the  shape  of  excess  absorption  will  be  given, 
following  consideration  of  the  observed  temperature  dependence. 

Typical  results  in  the  low  and  high  frequency  regions  are  given  in  Figure 
6.  The  excess  absorption  data  are  plotted  in  this  way  to  find  the  value  of 
the  energy  E  that  characterizes  the  temperature  dependence;  this  in  turn  can 
be  compared  with  the  estimated  dimer  energy  of  formation.  The  dotted  lines 
show  the  temperature  dependence  of  the  predicted  monomer  component  of  the 
absorption.  This  is  governed  mainly  by  the  Boltzmann  distribution  of  the 
energy  levels  from  which  nearby  absorption  transitions  arise;  the  slope 
therefore  switches  sign  from  one  side  of  the  monomer  pure  rotation  band  to 
the  other.  The  temperature  dependence  of  the  excess  absorption  is  more  nearly 
uniform.  Values  of  E  given  in  Figure  7  are  near  the  dimer  energy  of  formation 
estimated  from  the  temperature  dependence  of  the  second  virial  coefficient 
(Bolander  et  al.  1969;  Bohlander,  1979). 

5.  Discussion 


There  are  in  general  terms  two  possible  ways  of  interpreting  the  dis¬ 
crepancy  between  the  shape  of  observed  excess  absorption  and  that  of  predicted 
dimer  absorption: 
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(1)  Dimers  at  normal  temperatures  may  have  significantly  different  mean 
structures  than  they  have  at  low  temperatures,  or 

(2)  The  contributions  by  dimers  to  the  excess  absorption  may  be  less  Im¬ 
portant  than  possible  contributions  by  unbound  molecular  pairs,  l.e.t  there 
may  be  additional  collision  broadening  of  monomer  lines. 

Although  hydrogen-bonded  complexes  are  floppy  species.  It  was  not  entirely 
unreasonable  to  try  the  assumption  that  water  dimers  vibrate  harmonically  with 
small  amplitudes.  Thomas  (1975)  has  found  vibration  bands  of  the  complex  H^O — 
HF  which  were  easily  recognizable  despite  some  anharmoniclty.  Of  course,  ^0 — 
HF  has  a  much  higher  binding  energy  than  the  water  dimer  (an  energy  of  forma¬ 
tion  of  ~  0.3  eV  Instead  of  0.12  eV) ,  and  this  may  mean  that  barriers  to  in¬ 
ternal  rotation  are  also  higher.  If  water  dimers  have  low  barriers  to  internal 
rotation,  this  type  of  Internal  motion  could  be  more  important  than  vibration, 
and  the  mean  dimer  structure  could  be  significantly  different  than  it  is  at 
low  temperatures.  In  studies  of  the  more  weakly  bound  van  der  Vaals*  complexes, 
examples  have  been  found  of  vibrating  species,  of  internally  rotating  species, 
and  of  intermediate  cases  (Ewing  1976) .  According  to  the  theoretical  study 
made  of  barriers  of  internal  rotation  by  Owicki  et  al.  (1975),  mentioned  earlier, 
there  are  five  barriers  and  three  of  them  have  heights  less  than  100  cm  \ 

Since  1/2  kT  is  of  this  order,  it  is  likely  that  internal  rotation  is  important 
at  normal  temperatures.  Further  exploration  of  the  barriers  with  molecular 
orbital  methods  would  be  of  interest. 

The  other  principal  assumption  we  have  made  in  calculating  theoretical 
dimer  absorption  is  that  the  second  virial  coefficient  of  water  vapor  is 
entirely  due  to  the  presence  of  dimers,  and  on  this  assumption  we  have  esti¬ 
mated  the  concentration  of  dimers.  It  is  interesting  that  the  excess  absorp¬ 
tion  derived  from  observations  is  similar  in  magnitude  to  the  predicted  dimer 
absorption.  However,  the  principle  of  spectroscopic  stability  can  be  applied 
here  according  to  rules  given  by  Gordon  (1963),  and  one  finds  that  any  pair 
of  water  molecules,  whether  they  are  bound  or  unbound,  will  have  about  the 
same  integrated  absorption  cross-section  at  frequencies  in  the  far  infrared 
(below  the  frequencies  of  the  vibration  bands  of  the  monomer).  It  follows 
that  we  are  not  able  to  say  from  our  measurements  whether  the  second  virial 
coefficient  and  the  excess  absorption  are  due  to  bound  or  unbound  species. 
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Present  theories  discount  the  importance  of  unbound  pairs  at  normal 
temperatures,  but  there  are  some  simplifications  made.  For  molecules  with  a 
spherically  symmetric  form  of  potential  energy,  Stogryn  and  Hirschfelder 
< 1959)  have  shown  that  dimers  are  responsible  for  nearly  all  of  the  second 
virial  coefficient  when  the  dimer  binding  energy  exceeds  kT.  In  the  case 
of  water,  the  energy  of  dimer  formation  is  more  than  4kT,  but  it  is  not 
known  whether  the  generalization  of  Stogryn  and  Hirschfelder  can  be  extended 
to  the  case  of  molecules  of  low  symmetry,  such  as  water.  Calculations  of 
dimer  concentrations  may  also  be  done  from  a  knowledge  of  the  intermolecular 
potential  energy  (Viktorova  and  Zhevakin  1971;  Braun  and  Leidecker  1974;  Lane 
1975;  Lie  and  Clement!  1976:  Bohlander  1979).  While  reasonable  agreement 
with  experimental  values  of  the  second  virial  coefficient  have  been  obtained 
in  recent  calculations,  much  uncertainty  remains  about  the  reliability  of 
the  approximations  made.  Therefore,  the  possibility  of  a  significant 
contribution  by  unbound  pairs  to  excess  absorption  cannot  be  ruled  out. 

6.  Conclusions 


(1)  We  have  found  that  excess  absorption  by  water  vapor  is  similar  in 
magnitude  to  predicted  dimer  absorption  and  has  a  similar  temperature 
dependence. 

(2)  There  are  large  discrepancies  between  the  shape  of  excess  absorp¬ 
tion  and  that  predicted  for  dimers. 

(3)  Further  theoretical  work  is  needed  on  internal  rotation  in  dimers 
and  on  the  theoretical  concentrations  of  dimers  at  normal  temperatures. 
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Abs  tract 

The  near -millimeter  wavelength  region  (3.2mm  -  0.3mm)  is  being  investigated  for  military 
systems  applications  during  adverse  weather  and  in  the  presence  of  smokes,  dust  and  other 
particulate  clouds.  The  use  of  near-millimeter  wavelengths  (NMMW)  has  advantages  and  dis¬ 
advantages  relative  to  the  use  of  the  infrared  and  microwave  regions.  The  atmosphere  is  a 
dominant  factor  in  determining  the  operation  of  military  systems  in  the  NMMW  region. 

Systems  currently  under  consideration  for  NMMW  applications  include  beam  rider  and  command 
guidance,  missile  plume  detectors,  low-angle  tracking  radars,  terminal  homing  systems, 
target  acquisition  radars,  fuze  systems,  quasi -imaging  radars  and  hybrid  (1R/NMMW)  systems. 
Recent  NMMW  technological  developments  (e.  g.,  sources,  receivers,  components,  phenomenology 
and  measurements)  have  been  advancing  at  a  rapid  pace  to  meet  system  needs. 

I nt reduction 

Weapon  guidance  systems  operating  in  the  visible/infrared  spectral  regions  have  been 
successfully  developed  for  tactical  use  over  a  broad  spectrum  of  engagement  scenarios. 

These  include  direct  systems  employing  beam  riders,  command  guidance,  semi-active  and 
active  guidance  modes  which  have  demonstrated  effective  operation  in  clear  weather,  but 
ate  Inoperable  in  certain  severe  atmospheric  conditions  and  in  the  presence  of  smokes, 
dust  and  other  particulate  clouds.  Since  molecular  absorption  in  the  atmosphere  is  low 
over  broad  spectral  windows  throughout  the  visible/infrared,  the  primary  attenuation  mode 
is  Mie  scattering  by  water  droplets  in  haze  and  fog,  and  by  particulates  in  smokes  and 
other  aerosols. 

On  the  other  hand,  the  near-millimeter  wavelength  region  is  affected  less  by  adverse 
weather,  smokes,  and  aerosols  than  the  vtsible/3U  wavelength  regions.  Rain  attenuation 
of  NMM  wavelength  radiation  is  comparable  to  rain  attenuation  in  the  1R  region,  because 
or  its  Improved  transmission  under  adverse  environmental  conditions,  near-millimeter 
wavelengths  are  being  considered  for  military  applications.  The  NMMW  spectral  region  will 
not  prove  to  be  an  ideal  operational  region  for  all  systems  applications,  but  will,  in 
several  cases,  provide  a  compromise  adverse  weather  capability  for  limited  range  appli¬ 
cations  with  resolution  better  than  that  of  the  microwave  region. 

I  it  recent  years,  millimeter  and  near-millimeter  wave  systems  applications  have  been  the 
•ubtects  of  workshops',  reports-  and  study  panels^.  This  paper  will  draw  upon  material 
presented  in  these  references,  with  the  majority  of  Information  originating  from  the  Harry 
Diamond  Laboratories  study-'.  A  primary  purpose  of  many  or  these  investigations  has  been 
the  identification  of  NMMW  systems  technology  requirements,  some  of  which  will  be  briefly 
discussed  following  a  summary  of  potential  applications. 

Advantages  and  disadvantages  of  NMMW  legion 

The  near-millimeter  wavelength  region  has  usually  been  dismissed  in  the  past  as  a  region 
Of  high  attenuation,  relative  to  the  microwave  region,  and  as  a  region  of  poor  angular 
resolution  relative  to  t  tie  vislble/lR  wavelength  regions.  However,  several  factors  enter 
into  the  evaluation  of  the  wavelength  region  which  is  chosen  for  operation  of  a  particular 
system  during  adverse  weather.  ft  will  be  seen  in  t  lie  next  section  that  transmission  in 
adverse  weather  improves  considerably  in  going  from  visihle/TR  wavelengths  to  NMM 'MM 
wavelengths,  and  further  improves  at  centimeter  wavelengths.  On  the  basis  of  this  com- 
p.n  Ison,  it  would  appear  that  centimeter  or  longer  wavelengths  would  be  the  most  appropriate 
foi  systems  operation  in  adverse  weather.  Moreover,  it  lias  been  shown  that  transmission 
through  smokes  and  aerosols  also  follows  a  pattern  of  decreasing  with  increasing  wavelength. 
Aimpoint  accuracy  requirements,  however,  tend  to  eliminate  the  centimeter  wavelength  region 
on  the  basis  of  a  number  of  practical  considerations.  Tor  example,  the  required  accuracy 
for  target  tracking  or  target  designation  usually  tails  within  the  range  of  0.1  to  0.5  mi  ad 
resulting  in  antenna  dimensions  of  about  0.5  to  1  m  for  millimeter  wavelengths  and  5  to  10 
meters  for  centimeter  wavelengths.  This  consideration  alone  lilies  out  a  centimeter  wave¬ 
length  system  for  applications  where  the  weapon  system  platform  would  be  vehicle-  oi  heli¬ 
copter-mounted.  Analysis  of  NMMW  systems  requires  then  that  several  factors  be  weighed 
in  determining  their  potential  in  currently  conceived  tactical  applications.  Tables  1  -  4 
list  advantages  and  disadvantages  of  the  NMMW  spectral  region  relative  to  microwave  and 
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infrared  wavelengths.  Some  of  the  listed  advantages  might  be  questioned;  thus,  whereas  the 
best  of  RF  and  optical  techniques  can  be  combined  at  NMM  wavelengths,  tins  can  also  be 
considered  a  necessity  because  of  the  nature  ot  the  spectral  region.  On  the  other  hand, 
disadvantages  such  as  component  deficiencies  and  lack  of  data  will  be  reduced  ot  removed 
as  work  progresses  at  NMM  wavelengths. 

Whereas  several  of  the  NMMW  disadvantages  can  be  expected  to  diminish,  the  limiting 
factor  for  systems  applications  is  the  atmospheric  effect,  which  cannot  be  avoided.  The 
factor  e'"^.  where  a  is  the  atmospheric  attenuation  and  R  is  the  weapon- target  range,  is 
present  in  all  propagat ion  expressions.  In  the  analysis  and  discussions  of  NMMW  systems, 
their  performance  is  dominated  by  their  capability  to  operate  under  various  atmospheric 
conditions.  Since  the  malor  role  for  NMMW  tactical  systems  will  be  during  inclement 
weather  or  in  the  presence  of  smoke,  the  characteristics  of  a  NMMW  system,  when  such  con¬ 
ditions  prevail,  are  extremely  important. 

Table  1.  Advantages  for  NMMW  Region  Relative  to  Microwave  Region 

1.  Greater  Resolution 

2.  Smaller  Beam  Angle  t>  for  Given  Antenna  Diameter,  D,  and  Conversely  Smaller  D  for  Given 

"B  “ 

1,  Reduced  Multipath  Potentially  Improved  Low-Angle  Tracking 

4.  Low  Off-axis  Detectability  Providing  High  Security 

5.  Covertness  Due  to  Exponential  Eall-off  With  Range  in  High  Attenuation  Regions 

6.  Tracking  Through  Plumes  Improved  Over  Microwave  Systems 

7.  Smaller  Components  Allowing  More  Compact  On-board  Missile  Systems 

8.  Clutter  More  Diffuse,  Doppler  Shift  is  Greater,  Glint  Should  be  Less 

9.  Integration  of  Hybrid  IR/NMMW  Systems  Possible 

10.  Countermeasures  More  Difficult 

Table  2.  Advantages  for  NMMW  Region  Relative  to  Optical  Region 

1.  Improved  Transmission  in  Smoke  and  Fog,  Providing  Better  Low  Visibility  Operation 

2.  Harder  to  Jam;  More  Covert 

3.  Improved  Eye  Safety 

4.  Better  Coherent  Receiver  Techniques 

5.  Source  Stabilization  More  Easily  Performed 

6.  Potentially  Lower  Cost 

7.  Reduced  Background 

8.  Both  RF  and  Optical  Technology  Applicable 

Table  3.  Disadvantages  for  NMMW  Region  Relative  to  Microwave  Region 

1.  Poor  Heavy  Rain  Transmission 

2.  Atmospheric  (clear)  Absorption  is  Higher 

3.  Larger  Rain  Backscatter 

4.  Current  Receiver  Noise  Figures  are  Poorer 

5.  Higher  Precision  Manufacturing  Required 

6.  Solid  State  Source  Efficiencies  Fall  as  Frequency  Squared 

7.  Poorer  Source  Stability 

8.  Inferior  Waveguide  Power  Handling  Capability 

9.  Higher  Cost  for  NMMW  Systems 

Table  4,  Disadvantages  for  NMMW  Region  Relative  to  Optical  Region 

1.  Lower  Resolution 

2.  Higher  Glint 

3.  Lack  of  Relevant  Data 

4.  Video  Detection  Less  Sensitive 

5.  Component  Technology  Currently  Worse 

6.  NMMW  Laser  Sources  Less  Efficient 

7.  Solid  State  and  Tube  Sources  Chirped,  Less  Monochromatic  Than  Optical  Lasers 

8.  Components  Currently  Larger  and  Heavier 

Atmospheric  o ff ects 

Whereas  the  NMMW  region  offers  an  advantage  over  optical  systems  during  adverse  weather, 
clear  weather  NMMW  propagation  suffers  greater  attenuation  than  experienced  in  the  1R 
visible  region,  and  is,  in  addition,  for  most  systems  restricted  to  spectra l  windows.  The 
NMMW  region  has  stronq  Oj  and  11^0  absorption  lines  with  water  the  malor  absorber.  Between 
these  broad  linos  lie  the  transmission  windows.  Figure  1  shows  the  horizontal  attenuation 
across  the  NMMW  region  from  100  GHz  to  1000  GHz  at  sea  level  and  4  km  altitude'.  Also 
shown  is  the  weak  0j  absorption  at  sea  level.  Because  of  the  strong  absorption  at  the 
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shot  t  wavnlcngt  h  end  ol  tins  sped  ra  1  region,  t.i.-l  ir.tl  ground-to-ground  or  ground-to-ait 
apr  l  v  cut ions  .tii1  eont  ined  to  the  lonvii't  wavelength  windows  which  it  re  centered  About  the 
f  o  l  lowi  nq  wovelenut  Its  : 


Wavelength  (mm) 

Horizontal  Attenuation  ( dn/km) 

('em  t  h  At  t  enitat  ion  (dl>' 

3 

0.  I 

1.25 

2.  1 

0.5 

0.0  1 

1 .  3 

l  .75 

2.  1 

0 . 88 

9.0 

9.9 

0.72 

17 

22 

Included  in  this  listing  me  approximate  horizontal  attenuation  (dn/km)  for  5.>>1  q/m5 
of  HjO  at  sea  level  and  the  zenith  attenuation  (dB)  for  7.5  g/m'  of  HjO  at  sea  level.  A 
mot«'  comprehensive  indication  ot  .the  horizontal  attenuation  across  the  elect t omaunet ic 
spectrum  (3  cm  -  0.3  urn)  is  given  in  Figure  2.  A  clear  air  curve  for  HjO  density  of 
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Figure  1.  Spectral  plots  of  the  attenuation  by  the  (  1  o b 2 )  tt.^.  Standard  Atmosphere  at  s\'a 
level  and  4  km  altitude.  The  water  vapor  density  is  5.  91  u  m  at  sea  level  and  1.10  g.m’  at 
4  km  altitude.  The  lower  curve  represents  Os  only  at  sea  level.  For  comparison,  the  atten¬ 
uation  tor  this  same  model  atmosphere  is  approximately  0.2  dB/krn  in  the  10  micrometer  window 
and  less  than  0.1  uB/km  near  3.8  micromet ers[  4 ] . 

7.5  u  nt'  at  ,'0’C  shows  the  high  attenuation  in  the  submillimeter  region  surrounded  by  the 
low  microwave  and  IR/visiltle  attenuation.  From  the  fog  and  rain  curves  of  Figure  2,  the 
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Figure  2 .  Attenuation  by  atmospheric  gases,  rain  and  fog [5]. 

advant  age  of  NMMW  systems  during  heavy  fog  and  the  comparable  NMMW  and  IK  at  t  enuat ten  dm  tng 
rain  ate  demonstrated.  A  similar  figure  #  Figure  \  shown  increasing  attenuation  for  clouds 
with  decreasing  wavelength*  indicating  an  advantage  for  satellite  or  at tcra! t - to-gr 'ound  s\  *- 
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terns  applications.  Near  millimeter  wave  attenuation  by  aerosols  is  highly  dependent  on  the 
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Figure  3.  Summary  of  sea-level  atmospheric  attenuation[6] . 

type  of  aerosol  and  its  particle  size  distribution  and  concentration.  Table  5  shows  the 
size  range  and  concentration  of  water  droplets  in  rain,  fog  and  clouds.  It  can  be  seen  from 
this  table  that  the  larger  size  range  of  rain  droplets  (7-700  urn)  would  cause  significant 
scattering  losses  at  near  millimeter  wavelengths  were  its  concentration  not  extremely  low 
compared  to  the  concentration  of  water  droplets  in  fogs  and  clouds.  For  both  rain  and  fog, 
the  attenuation  is  an  almost  linear  function  of  bulk  water  density  with  the  frequency  de¬ 
pendence  for  attenuation  in  rain  being  negligible  over  the  NMM  spectral  region  from  94  to 
340  GHz. 


Table  5.  Comparative  Water  Droplet  Size  and  Density 
_ For  Rain,  Fog,  and  Clouds  [8,  9,  101 _ 


Sin*  Range  (pm) 
Droplet  Density  (m  3) 
Water  Density  (g  m3) 


Light 

*  1  mm  'hr) 


7  -  100 
350 
0.04 


Moderate 
( r  *  4  mm  tor) 


10  -  300 
500 
0.17 


Heavy 

(r  «  25  mm  tor) 


10  -  500 
700 
1.0 


Cloudburst 
(r  *  100  mm  tor) 


50  -  700 
1,250 
4.2 


Fog 


Siic  Range  (pm) 
Droplet  Density  (m  3) 
W  ater  Density  (g  m3) 


Thin  (vis.  -  300  m) 


2  *  10‘ 
0.023 


Thick  (vis.  -  50  m) 


0.02  -  0.2 
.3  v  1011 
0.35 


Si/c  Range  (pm) 
Droplet  Density  (m  3) 
W  ater  Density  (g  m  ) 


Cumulonimbus 

Cumulus  Conge  stus 

Fair  Weather  Cumulus 

Strato  Cumulus 

6-12 

3-5 

- 

- 

1.6  *  103 

6.0 

2.5 

0.5 

m 

For  aerosols  resulting  from  smokes,  dust  or  other  battlef ield-generatod  debris,  the 
attenuation  is  primarily  due  to  scattering,  and  secondarily  to  absorption.  Preliminary 
indications7  are  that  propagation  through  existing  screening  smokes  is  very  high,  and  that 
the  development  of  effective  smokes  for  the  NMMW  reqion  is  questionable  because  of  the 
difficulties  involved  with  generating  and  dispensing  smokes  with  the  appropriate  particle 
size  distribution  and  settling  rate.  Attenuation  of  millimeter  radiation  by  ground  ex- 


plosions  of  munitions  has  been  shown  to  be  relatively  high  for  short  periods  of  time7. 
Although  the  attenuation  cleared  much  faster  for  NMM  waves  than  for  the  IR/optical  bands, 
nevertheless,  loss  of  NMMW  track  or  acquisition  of  a  target  could  occur  during  a  barrage 
or  critically  placed  explosion  so  that  these  effects  must  be  considered  in  greater  detail 
in  the  future. 


It  is  important  at  this  point  to  demonstrate  millimeter  wave  propagation  under  condi¬ 
tions  of  high  absolute  humidity  and  high  bulk  water  content,  since  these  are  the  condi¬ 
tions  under  which  all  military  systems  are  expected  to  operate.  Examples  are  given  for 
the  three  cases:  clear  weather  (high  absolute  humidity),  rain  and  fog. 

Clear  weather 


In  high  visibility  conditions  with  a  lack  of  precipitation  and  fog,  the  attenuation  of 
near  millimeter  radiation  in  dB/km  is  directly  proportional  to  absolute  humidity.  This 
linear  relationship  is  shown  in  Figure  4  on  a  plot  of  temperature-humidity  data  which  was 
abstracted  from  a  psychrometric  chart.  Table  6  shows  the  two-way  path  loss  for  a  range 
of  3  km  at  frequencies  of  94,  140  and  220  GHz  for  various  atmospheric  conditions.  Note 
that  the  attenuation  of  220  GHz  radiation  is  very  high  for  warm,  humid  conditions  although 
the  total  path  loss  is  only  18  dB  (3  dB/km)  for  the  standard  atmosphere  (T  =  300'K, 
o  »  7.5  g/m^ ) . 


STANDARD  ATMOSfU 
T  ■  J00°K..i  •  7  5 f /m" 


Figure  4.  Atmospheric  absorption  by  water  vapor 


Table  6.  Two-Way  Path  Loss  In  dB  For  Atmospheric  Attenuation  Over  A  Range  of  3  km 


r .  wk 

T  •  WF 

r  *  40* 

F 

30', 

50\ 

30'1 

:>o'\ 

so'l 

30H 

r,o«; 

S0‘* 

3. 1 

3.2 

S,  4 

1.4 

1.9 

3.0 

0.5 

o.o 

1.4 

13.2 

21.0 

35.4 

5.1 

0.3 

15.0 

2.7 

3.3 

G.2 

jG.  4 

13.  * 

. 

70,  Jt 

11.4 

1*.C 

30.0 

4.5 

G.  3 

12.3 

The  impact  of  these  attenuation  figures  on  system  design  can  be  shown  by  solving  the 
radar  equation  for  a  typical  set  of  system  parameters;  the  radar  range  equation  is  of 
the  form 


S/N 


pptG?*2°T 
(  4  n  )  3R‘*KTNfLT 


(1) 


where  Pp  »  peak  power  of  the  transmitter;  x  =  pulse  width  of  the  transmitter  (it 

has  been  assumed  that  the  receiver  bandwidth  =  1/t) ;  G  *  gain  of  the  trans¬ 
mitting/receiving  antenna;  1  =  free  space  wavelength;  of  =  radar  cross  sec¬ 
tion  of  the  target;  R  ■  range  to  target;  K  =  Boltzmann's  constant;  Nf  =  noise 
figure  of  the  receiver,  and  =  total  radar  system  loss  (atmospheric  &  signal 
processing  +  waveguide  and  components) . 

Assuming  the  following  parameters: 


t  =  200  nsec.;  \  =  1.3  mm  (230  GHz);  G  =  47.6  dB  (6"  aperture  with  n  »  0.5); 
jT  »  50  m  ,•  KT  =  4.14  x  10~21  (t  =  300”K);  Ls  -  signal  processing  loss  «  3  dB; 
Ly  «  waveguide  and  component  loss  =  8  dB;  =  12  dB. 


we  may  write 


S/N (dB)  *  PT(dBW)  -40  log  R(km)  -  2n (dB/km) R(km)  +  15.2 


1 


Figures  5  and  6  show  S/N  vs  R  for  different  values  of  PT  for  atmospheric  conditions  of  T  »  27°C  and  p  =  7.5 
g/m3  (standard  atmosphere)  and  T  =  30°C  and  p  =  28  g/m3  (relative  humidity  *  80%).  For  these  two  cases, 
Table  7  shows  the  peak  power  versus  range  required  for  S/N  =  14  dB  (P,  =  6.9,  Pf  »  10-7) .  Note  that  for 
large  a  ,  the  S^l  vs  R  curves  are  very  steep,  and  only  small  changes  in  range  occSr  for  large  changes  in  By 
at  a  given  S/N.  With  a  maximum  peak  power  of  10  kW,  for  example,  the  high  visibility  range  will  be  limited 
to  about  1.5  km  on  a  warm  humid  day,  and  to  about  3  km  on  a  warm  dry  day.  These  calculations  are  for  single 
pulse  transmission  only. 
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Figure  6.  Performance  of  a  230  GHz  pulsed  system  on  a  warm,  humid  day. 


Table  7,  Required  Power  To  Achieve  S/N  =  14  dB 


Warm,  Dry  , 

Warm,  Humid 

T  =  27°C,  p  =  7.5  g/m 

T  «  30  °C ,  p  ■  28 

Range  (km) 

Pt(KW) 

Pt(KW) 

1.0 

0.1 

1.25 

1.0 

1.5 

7.2 

1.75 

47.0 

2.0 

0.2 

284.5 

2.25 

0.4 

2.5 

0.9 

2.75 

1.8 

3.0 

_  3.6 

Although  the  performance  of  a  system  operating  in  the  230  GHz  window  degrades  rapidly 
in  warm  humid  weather,  in  many  parts  of  the  world  these  conditions  are  present  for  only 
a  small  percentage  of  the  time.  An  example  is  the  climatological  conditions  of  Western 
Europe.  Figure  7  shows  the  attenuation  for  three  near-millimeter  frequencies  based  on 
the  mean  maximum  daily  temperature  and  relative  humidity  for  Bayreuth,  West  Germany. 
These  data  predict  that  the  mean  attenuation  during  the  summer  months  will  be  less  than 
7  dB/km,  and  for  this  figure,  a  range  of  2.0  km  would  be  possible  with  a  10  kw  pulse 
source. 


_  I  i m f  1  ----- . 
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Figure  7.  Millimeter  wave  atmospheric  absorption  by  water  vapor  -  Bayreuth,  West  Germany. 

Rain.  Although  differences  exist  in  predicted  rain  attenuation  for  the  NMMW  region1*,  most 
models  predict  the  following  relationships: 

a(dB/km)  «  ArB,  (2) 

where  A  and  B  are  constants  for  the  frequency  region  94  -  340  GHz,  and  r  is  the  rain  rate 
in  millimeters  per  hour.  Based  on  the  available  data,  the  relation 

a  (dB/km)  «  r0-85  (3) 

provides  a  reasonably  close  fit  throughout  the  NMMW  spectrum.  As  might  be  expected,  a  warm 

rainy  day  will  be  a  problem  for  a  system  at  220  GHz  because  of  strong  absorption  by  both 

water  vapor  and  bulk  water.  For  example  at  T  *  26.7*C  (80'F)  ,  and  r  «  4  mm/hr.  (a  moderate 
rain  rate),  the  attenuation  is  11  dB/km  (water  vapor)  plus  3  dB/km  (rain)  or  14  dB/km  total. 
Figure  8  shows  the  S/N  versus  range  at  three  power  levels,  and  again  it  is  evident  that  for 
a  given  value  of  S/N,  large  increases  in  power  result  in  small  increases  in  range.  For 

this  atmospheric  condition,  a  peak  power  level  of  about  50  kW  would  be  required  for  a  range 

of  1.5  km  at  S/N  »  14  dB  whereas  a  range  of  about  1.2  km  would  be  possible  at  the  same 
value  of  S/N  with  only  5  kW  peak  power.  An  order  of  magnitude  increase  in  power  only  re¬ 
sults  in  about  0.3  km  increase  in  range  for  this  condition. 

Fog .  Attenuation  of  NMMW  radiation  by  fog  is  also  dominated  by  bulk  water  absorption, 
although  a  slight  frequency  dependence  is  also  present.  A  good  fit  to  available  fog 
attenuation  data  is  provided  by  the  expression 


a  (dB/km) 


0.035  plg/m3)  f1-03(GHz) 


This  relationship  along  with  the  optical  visibility  is  shown  in  Figure  9.  For  a  30  meter 
visibility  radiation  fog,  with  liquid  water  content  of  0.75  g/m-’,  the  attenuation  by  water 
absorption  at  230  GHz  is  approximately  9  dB/km.  As  shown  in  Figure  7,  the  mean  maximum 
daily  temperature  of  Bayreuth,  West  Germany  during  the  winter  months  (the  season  ol 
greatest  incidence  of  fog)  is  about  0-4*C;  therefore  a  reasonable  fog  scenario  is  30-100  m 
visibility  with  temperature  0-4 *C,  for  which  the  attenuation  at  230  GHz  would  be  approxi¬ 
mately  4  dB/km  to  11.5  dB/km.  These  conditions  have  been  bracketed  by  the  calculations 
presented  in  Figures  5  and  8.  It  is  significant  to  point  out  that  a  100  m  visibility  fog 
at  T  »  4*C  (40*F)  results  in  an  attenuation  of  only  4.5  dB/km  at  230  GHz  compared  with 
about  50  dB/km  in  the  8-12  urn  band  and  about  150  dB/km  at  a  wavelength  of  0.6  urn. 
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where  T„  is  the  temperature  of  the  earth- or  qround-based  target  and  T(z)is  the  temperature 
of  a  stratum  of  atmosphere  of  thickness  dz  located  at  altitude  Z.  The  terms  of  the  form 
t(Z,,  are  the  optical  depths  between  altitudes  Z,  and  Z,  at  angle  0,  and  R  is  the 

reflectivity  of  the  earth  or  target.  The  optical  depths  tare  given  by 


VISIBILITY  (m) 


Table  8.  Summary  Of  Propagation  Data 


Attenuation,  a  (dB/km) 
Water  Density,  p  (g/m’) 


i  r-  o 


a  Clear 
Rel.  Humidity 


a  Fog 


a  Fog 


a  Cloud 


X  (pm) 

-  100Z 

T  =■  32"  F 
p  -  4.8 

68*F 

17.4 

Radiative 
Rv  -  400ra 
p  -  0.014 

Fog 

200 

0.038 

100 

0.11 

30 

0.71 

Advectlon  Fog 

400  200  100 

0.063  0.18  0.4 

a 

1 

Rain 

mm/hr 

4 

10 

Fair 

Weather 

Cum. 

Nimbo- 

Strat. 

1 

80 

500 

95 

570 

4 

120 

640 

10.6 

0.3 

1.2 

7 

20 

58 

373 

17 

63 

140 

1 

2.6 

6 

50 

500 

337 

50 

185 

0.6 

1.5 

4.3 

28 

2.5 

7.1 

15.8 

1 

3 

7 

3 

20 

72  It 

10 

37 

0.3. 

0.9 

2.6 

17 

1.0 

4.3 

9.6 

1 

3 

7 

2 

7 

880 

7 

24 

0.3 

0.7 

2.0 

14 

1.2 

3.5 

7.9 

1 

3 

7 

1.5 

6 

1300 

2 

6 

0.2 

0.5 

1.4 

9 

0.8 

2.3 

5.1 

1 

3 

7 

0.8 

4 

2300 

1 

3 

0.1 

0.2 

0.6 

4 

0.4 

1.0 

2.2 

1 

3 

8 

2 

3200 

0.2 

0.9 

0.1 

0.2 

0.5 

3.2 

0.3 

0.8 

2 

1 

3 

8 

1.5 

Notes: 

(!)  f°r°CLEAR 

at  other  Rel 

.  Hum. , 

scale 

down 

from  100Z 

given 

(2)  for  a 

fog  Situation,  T0TAL  - « 

CLEAR 

(RH  - 

■  100Z)  + 

c^,  (likewise 

for  clouds) 

(3)  for  a 

rain 

situation,  <*, 

TOTAL  “ 

“clear 

-  100Z)  + 

a 

R 

t(Z-^,  Z  2 , 


a(Z) 


sec  edZ 


(6) 


so  that  the  first  term  of  Equation  (5)  is  the  contribution  of  the  intervening  atmosphere 
between  the  radiometer  and  target,  the  second  term  represents  atmospheric  emission  re¬ 
flected  frcm  the  targetto  the  radiometer,  and  the  third  term  is  the  target's  emission  atten¬ 
uated  by  the  atmosphere  between  the  radiometer  and  target.  The  target  or  earth  background 
reflectivity  is  given  by  R,  the  target  (earth)  temperature  is  T£  and  the  atmospheric  atten¬ 
uation  is  a ( Z , v ) .  It  is  seen  from  these  relations  that  not  only  atmospheric  absorption 
by  02  and  H20  is  important,  but  hydrometeorite  absorption  and  scattering  and  atmospheric 
thermal  emission  are  also  determining  factors  for  the  radiometric  antenna  temperature, 
which  is  reduced  further  in  value  when  the  target  fill  factor  for  the  antenna  is  less  than 
one. 


Preissner^  has  shown  the  brightness  temperature  contrast  for  different  weather  con¬ 
ditions  and  different  materials  for  the  microwave  through  NMMW  spectral  region.  Figure  10 
demonstrates  the  brightness  temperature  contrast  between  vegetation  and  concrete  for 
three  different  altitudes  for  a  clear  standard  atmosphere  (p  _  =  7.5  g/m  ,  T  =  20°  C  and 
P  =  760  mmHg  at  sea  level) .  At  all  altitudes  given,  a  detectable  brightness  temperature 
contrast  exists,  except  possibly  at  230  GHz  for  3  km  and  8  km,  where  the  temperature  con¬ 
trast  is  on  the  order  of  3-5^  K  and  the  minimum  detectable  temperature  is  approximately 
0.5-3°  K  for  a  1  GHz  bandwidth  and  a  10  msec  time  constant.  Figure  11  shows  the  bright¬ 
ness  temperature  contrast  (concrete  and  vegetation)  for  a  radiometer  at  3  km  altitude  for 
several  atmospheric  conditions.  Rain  severely  degrades  NMMW  radiometric  response  and  fog 
of  visibility  of  100  m  or  less  results  in  considerable  loss  of  contrast.  A  more  signif¬ 
icant  indication  of  radiometric  capability  to  detect  a  military  target  is  the  brightness 
temperature  contrast  for  metal  relative  to  vegetation.  Figures  12a  and  12b  show  contrasts 
for  metal,  water  and  concrete  relative  to  vegetation.  For  metal,  detectable  contrasts 
are  observed  for  all  conditions  except  the  25  mm/hr  rain.  It  must  be  emphasized,  however, 
that,  for  metal,  reflectivity  of  1  has  been  used  and  no  fill  factor  reduction  has  been 
employed.  Both  conditions  are  unlikely.  Systems  applications  require  further  consider- 
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ation  of  these  effects  with  inclusion  of  antenna  pattern  and  fill  factor  effects 


Very  little  data  exist  for  atmospheric  turbulence  in  the  NMMW  region;  fluctuation  of 
intensity  and  signal  angle  cf  arrival  for  NMMW  systems  have  been  studied  by  Armand  et 
al.  and  by  Snider,  Wiltse  and  McMillan^1*.  Based  on  the  approach  first  proposed  by 
Armand,  the  peak-to-peak  intensity  fluctuations  expected  for  a  210  GHz  system  under  near 
standard  conditions  were  calculated  to  be  1.16  dB,  and  the  corresponding  angle-of -arrival 
fluctuations  were  determined  to  be  0.1  mrad.  in  practice,  observed  intensity  fluctuations 
at  lower  frequencies  have  occasionally  been  observed  to  be  much  larger,  and  little  data 
are  available  for  the  angular  fluctuations.  Much  work  needs  to  be  done  before  the  effects 
of  turbulence  on  NMMW  systems  have  understood,  but  the  meager  results  available  to  date 
indicate  that  it  gives  system  degradation  of  a  magnitude  which  must  be  considered. 


The  approximations  discussed  thus  far  demonstrate  some  of  the  atmospheric  effects  to 
be  considered  in  systems  trade-offs  between  the  NMMW  and  IR./visible  regions  but  do,  in 
no  way,  describe  the  complexity  that  could  occur  at  NMM  wavelengths.  Thus,  when  the 
limited  experimental  data  for  HjO  are  compared  with  existing  theories,  the  measured  values 
in  the  window  regions  invariably  exceed  theoretical  values.  When  measured  values  a>. e  com¬ 
pared  with  the  monomer  spectrum  calculated  with  the  Van  Vleok-Weisskopt ‘  '  or  Grossly  line- 
shape,  the  inadequacy  of  the  theory  is  evident  from  Figure  11  in  which  the  dashed  curve  B 
represents  an  empirical  "continuum"  absorption  which  must  be  added  to  bring  the  theoretical 
results  into  coincidence  with  the  experimental  results'.  The  empirical  correction  term 
has  been  given  in  the  format 


where  p  is  the  water  density  (g/m3) ,  T  is  the  atmospheric  temperature,  °K,  F  is  the  pres¬ 
sure  and  v the  frequency.  The  discrepancy  between  measured  and  calculated  water-vapor 
absorption  closely  follows  this  empirical  correction  term  throughout  the  NMMW  region. 
Various  causes  have  been  postulated  for  this  excess  or  anomalous  absorption,  but,  thus  far 
the  source  is  not  understood.  Great  anomalies  have  been  observed  in  high  humidity,  cold 
conditions  and  in  fog1  .  Further  experimentation  is  needed  in  this  area  of  atmospheric 
effects . 
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Figure  11.  Spectral  plots  of  the  attenuation  by  atmospheric  H ,0  at  sea  level.  Curve  A 
represents  a  combination  of  theoretical  and  experimental  results  for  an  H  ,0 
density  of  5.61  g/m-*.  At  a  fixed  temperature  and  1  atm  total  pressure,  the 
attenuation  is  approximately  proportional  to  the  H ,0  density.  Curve  R  cor¬ 
responds  to  an  empirical  continuum  that  is  added  tii  theoretical  results  to 
provide  agreement  with  the  experimental  results.  The  transrai t tance  scale  on 
the  right-hand  side  corresponds  to  a  1  km  path[4]. 

Cons idorab le  space  has  been  devoted  in  this  paper  to  a  discussion  of  at mospher ic 
effects  on  propagation.  However,  it  is  just  these  effects  which  are  the  most  important 
factors  in  determining  the  applicability  of  a  particular  military  system  in  the  NMM  wave 
length  region. 
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Potential  NMMW systems  app 1 i cations 


Several  NMMW  military  applications  have  been  investigated  in  recent  years.  Of  these 
some  have  not  shown  advantages  over  the  equivalent  system  operating  in  another  spectral 
region,  whereas  some  have  sufficient  promise  to  warrant  initiation  of  experimental 
development  of  prototype  systems.  Among  the  potential  systems  applications  which  have 
been  studied  are  the  following*'^ 


1.  Beamrider  systems 

2.  Terminal  homing 

3.  Target  designation  and  semi-active  homing 

4.  Command  guidance 

5.  Target  surveillance  and  acquisition 

6.  Active  NMMW  quasi- imaging 

7.  Low  angle  tracking 

8.  IFF  systems 

9.  Airborne  passive  imaging 

10.  Aircraft  detection  from  satellites 

11.  Boost  phase  plume  detection 

12.  Re-entry  applications 

13.  Secure  communications 

14.  Mine  detection 

15.  Obstacle  and  terrain  avoidance 

16.  Space  object  identification  (SOI) 

17.  Fuzing 

18.  Hybrid  (IR/NMMW)  systems 

The  details  of  the  investigation  of  most  of  these  systems  can  be  found  in  the  references 
and  only  a  selected  system,  the  NMMW  beamrider,  and  brief  statements  on  others  will  be  dis¬ 
cussed  in  this  section. 


Beamrider 


A  beamrider  guidance  system  '  '  is  defined  as  a  technique  for  guiding  missiles 

which  utilizes  a  beam  directed  into  space,  such  that  the  beam  axis  forms  a  line  along  which 
it  is  desired  to  direct  a  missile.  The  missile  contains  equipment  that  can  sense  the 
direction  and  magnitude  of  the  error  when  its  path  has  deviated  from  the  center  of  the  beam 
and  that  can  generate  guidance  error  signals  which  cause  the  missile  to  return  toward  the 
center  of  the  beam.  In  principle,  a  beamrider  system  can  be  employed  in  surface-to-surface 
surface-to-air,  air-to-air,  and  air-to-surface  roles,  although  it  is  potentially  more 
effective  against  slowly  moving  targets.  The  basic  elements  of  a  beamrider  system  are  the 
same  for  each  of  the  roles  but  the  stringency  of  the  requirements  placed  on  these  elements 
will  be  determined  by  the  particular  applications.  For  the  purposes  of  interest  here,  a 
surface-to-surface  anti-armor  role  has  been  emphasized  as  an  example  of  a  credible  appli¬ 
cation  of  NMMW  guidance. 


In  the  beamrider  concept  adopted  as  an  example  for  this  study,  the  target  is  acquired 
and  tracked  by  a  precision  near  millimeter  wave  (NMMW)  radar  with  a  conical  scan  antenna. 
The  missile  is  launched  toward  the  target  along  the  conical  scanning  tracking  beam,  which 
is  coded  to  provide  information  on  the  position  of  the  missile  relative  to  the  scan  axis 
which  defines  the  line-of-sight  (LOS)  to  the  target.  A  receiver  in  the  rear  of  the  missile 
detects  the  scanning  beam,  decodes  it  to  determine  its  position  with  respect  to  the  scan 
axis,  and  generates  the  appropriate  error  signals  to  keep  the  missile  in  the  center  of  the 
beam.  In  the  conical  scan  mode,  a  relatively  simple  PRF  coding  scheme  can  be  used  to 
develop  the  required  guidance  error  signals. 


Figure  14  shows  the  NMMW  beamrider  guidance  concept.  The  first  and  second  beams  are  the 
capture  beams  which  establish  line-of-sight  (LOS)  guidance  as  quickly  as  possible.  The 
basic  functions  of  the  beamrider  system  are  target  acquisition,  target  track,  missile 
capture  and  missile  guidance.  /■=--■ v_-  , 
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Figure  14.  NMMW  beamrider  concept 
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The  basic  system  configuration  considered  to  be  most  desirable  for  a  ground-to-ground, 
anti-armor  beamrider  system  requires  that  the  entire  system  be  a  self-contained,  crew- 
served  and  vehicle-mounted  package  with  an  antenna  aperture  diameter  less  than  one  meter. 
The  candidate  missile  would  be  of  the  TOW  or  SHILLELAGH  generic  class,  and  would  be  tube- 
launched  from  the  tracking/guidance  platform.  This  concept  would  require  an  integrated 
target  acquisition  capability  for  performing  target  handover  from  a  wide-area  battlefield 
surveillance  system,  and  for  establishing  an  autonomous  operating  capability  over  a 
limited  battlefield  sector.  The  acquisition  capability  would  of  necessity  require  the 
development  of  suitable  target  recognition  criteria  based  on  the  near  millimeter  wave 
target  signature. 


Nominal  system  parameters  relating  to  these  requirements  have  been  given  as  the 
following  : 


Operating  Range 
Capture  Range 
Weather  and  Environ¬ 
mental  Extremes 


Acquisition 
Missile  Data  Rate 


-  Maximum:  >  2  km;  Minimum:  <  0.5  km 

-  Maximum:  100  m 

-  Rain:  4  mm/hr;  Fog:  100  m  vis. 

Temperature:  90  -  100°F; 

Relative  Humidity:  80  -  100% 

Smoke:  Battlefield-generated  dust  and 

aerosols;  tactical  screening  agents 

-  Search  Sector:  5°  x  25°,  movable  scan  sector; 
Search  Time:  <  10  sec. 

-  50  Hz 


In  addition  to  these  characteristics,  a  NMMW  system  should  also  satisfy  the  general  re¬ 
quirements  of  high  reliability  and  maintainability,  compatibility  with  existing  equipment, 
capability  for  CM  hardening,  minimum  operator  interaction  and  training,  rapid  fire  capa¬ 
bility,  and  high  first  round  probability  of  kill. 


Based  on  these  requirements,  system  concepts  and  estimated  performance  for  a  NMMW  anti¬ 
armor  beamrider  system  have  been  addressed^.  However,  the  paucity  of  essential  data  for 
target  and  terrain  reflectivity  and  atmospheric  propagation  may  result  in  an  altering  of 
the  performance  calculations  when  those  data  become  available. 


The  beamrider  concept  has  been  successfully  implemented  in  the  infrared,  demonstrating 
good  clear  weather  performance.  Currently,  there  is  an  enthusiastic  interest  in  the  beam¬ 
rider  for  the  anti-armor  application.  The  U.  S.  Army  Ballistic  Research  Laboratory  (BRL) 
has  successfully  carried  out  a  feasibility  demonstration  o'?  beamrider  tracking  and  guidance 
at  140  GHz 20,  and  plans  to  repeat  these  experiments  at  217  GHz.  Efforts  have  also  been 
initiated  at  U.  S.  Army  MIRADCOM  to  define  and  verify  a  baseline  beamrider  missile  system 
for  anti-armor  applications. 


The  critical  performance  parameters  of  any  beamrider  system  are  related  to  (1)  high 
accuracy  centroid  tracking  of  the  target,  (2)  a  temporal  and  spatial  beam  structure  such 
that  the  missile  guidance  is  not  biased  by  the  terrain  effects,  and  (3)  a  method  of 
capturing  the  missile  at  a  range  close  enough  to  insure  adequate  damping  of  the  missile's 
angular  deviations  from  LOS  at  the  closest  tactical  range  of  interest. 


The  application  of  beamrider  technology  is  most  promising  in  the  guided  direct  fire 
anti-tank  role  with  the  radar  and  launcher  borne  by  a  land  vehicle.  The  system  as  pre¬ 
sently  envisioned  requires  a  relatively  large  antenna  0.6  m  diameter),  a  precision 
tracking  mount,  and  a  dual  beamwidth  antenna  configuration  to  accomplish  missile  capture. 

It  is  not  feasible  to  consider  a  helicopter-borne  weapon  system  at  this  time  because  of 
the  size  of  the  antenna/mount  configuration;  however,  it  may  be  possible  to  exploit  optical 
scanning  techniques  which  involve  stabilization  of  scan  components  to  obviate  the  require¬ 
ment  for  a  large  precision  tracking  mount.  The  anti-helicopter  role  could  be  considered 
if  the  target  exposure  time  is  adequate  to  acquire  and  to  accommodate  the  missile  time-of- 
flight.  In  the  anti-aircraft  role  the  end  game  maneuverability  requirements  on  the  missile 
would  be  severe  for  crossing  targets  or  for  targets  performing  high  G  maneuvers;  however, 
such  a  system  appears  to  be  feasible. 

In  considering  why  a  beamrider  system  should  be  investigated  as  the  initial  concept  for 
NMMW  guidance  feasibility  rather  than  another  scheme,  the  following  factors  are  influential: 

(1)  The  guidance  link  is  one-way,  so  that  relatively  simple  video  receivers  are  poten¬ 
tially  adequate  for  the  missile  link. 

(2)  Beamrider  guidance  has  been  demonstrated  at  IR  and  optical  wavelengths  in  opera¬ 
tional  situations  so  that  it  is  not  necessary  to  go  through  a  concept  demonstration  phase. 

(3)  Target  cross  section  is  high  because  the  tracker  looks  at  backscatter  rather  than 
at  off-axis  scatter. 
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(4)  Since  only  the  low  cost  missile  receiver  is  expended  in  firing,  the  cost  per  round 
could  be  significantly  less  than  that  for  a  missile  employing  an  active  seeker. 

(5)  The  missile  is  difficult  to  jam  since  the  receiver  looks  back  at  friendly  territory. 

(6)  It  is  possible  to  implement  a  multiple  wavelength  system  so  that  wavelength  opti¬ 
mization  is  possible  for  the  acquisition,  tracking  and  guidance  functions. 

(?)  The  airframe  could  be  optimized  aerodynamical ly  since  there  is  no  seeker  to  compli¬ 
cate  the  warhead.  Complexity  of  the  on-board  guidance  and  control  is  minimal. 

The  NMMW  beamrider  system  must  perform  several  functions  in  sequence,  with  the  capabi¬ 
lity  for  continuity  or  rapid  handover  from  one  function  to  the  next.  This  system  must  be 
capable  of  searching  for  and  acquiring  a  target  over  a  limited  f ield-of-view  (FGV) , 
followed  by  the  target  tracking  and  missile  guidance.  Handover  from  the  tracking  to  gui¬ 
dance  modes  must  include  a  missile  capture  phase  during  which  the  missile  is  launched 
toward  the  target  and  LOS  guidance  along  the  tracking  radar  beam  is  established.  For 
simplicity,  it  is  desirable  that  one  system  at  a  single  wavelength  perform  the  entire 
beamrider  operation;  however,  because  of  the  complexity  of  these  multiple  operations 
(target  acquisition,  tracking,  missile  capture  and  guidance) ,  the  optimum  configuration 
might  require  a  multi-wavelength  system. 

Several  approaches  for  implementation  of  the  beamrider  system  have  been  investigated2^. 
Performance  in  inclement  weather  is  the  most  important  consideration  for  a  NMMW  beamrider. 
For  94  GHz,  140  GHz  and  220  GHz,  an  example  of  the  transmitter  power  requirement  for 
target  acquisition  as  a  function  of  range  for  rain  rates  of  4  and  16  mm/hr  on  a  warm  day 
is  given  in  Figure  15.  The  following  parameters  are  used  in  the  calculations  for  source 
power : 


1  (m) 

G„  (dia.  ■  0.6  m;  n 
cl  (m2) 

S% 

kTB  (B  »  107  Hz) 

Nf  (dB> 


94  GHz 


140  GHz 


3.19  X  10"3 
0.5)  1.78  x  105 

30 
0.5 

3.1  (4.9  dB) 
4.14  x  10'14 
4 


2.14  x  10'3 
3.88  x  105 

30 

0.40 

6.2  (7.9  dB) 

4.14  x  10'14 
7 


220  GHz 

1.36  x  10*3 
9.60  x  105 
30 

0.16 

38.9  (15.9  dB) 
4.14  x  10-14 
15 


Figure  15.  Power  required  for  target  acquisition  on  a  warm  day  with  moderate  to  heavy  rain 
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Trade-offs  between  search  time,  PRF ,  and  antenna  beamwidth  must  Vie  considered  in  de¬ 
signing  a  beamrider  system.  If  a  separate  acquisition  system  is  employed,  the  required 
S/N  for  tracking  will  be  considerably  improved  by  the  integration  of  additional  pulses 
during  the  tracking  interval  which  is  long  with  respect  to  the  acquisition  target  dwell 
time. 

The  curves  of  Figure  15  show  that  only  the  44  GHz  system  is  capable  of  a  range  of  3  km 
with  a  power  of  1  kW  in  moderate  rainfall.  The  transmitter  power  requirements  lor  missile 
guidance  are  not  as  severe  as  those  required  for  acquisition  and  tracking  because  the 
guidance  link  is  a  one-way  path.  It  has  been  shown**  that  power  levels  required  for  acqui¬ 
sition  are  more  than  adequate  for  guidance. 

For  a  five  year  projection  of  system  performance,  it  would  seem  reasonable  to  anticipate 
improved  noise  figures  and  decreased  RF  losses  at  140  and  220  GHz.  These  projections  of 
improved  performance  are  listed  below: 


Frequency  (GHz) 

State 

-of-the-art 

5  Year  Projec 

t  UM1 

*‘RF 

NEP  N 

f 

lrf 

NEP 

Nf 

140 

4  dR  1. 

6  x  10‘12  7 

dB 

3  dB 

IQ'12 

6  dB 

220 

8  dB 

5  x  10'12  15 

dB 

3  dB 

io‘12 

6  dB 

reason  for  projecting 

the  same  paramo 

ters  for  both 

frequencies  is 

that  qua 

S 1 -opt 

components  should  show  very  little  frequency  dependence  over  the  range,  and  the  cutoff 
frequency  of  Schottky  barrier  mixers  is  presently  about  .1000  GHz. 

Target  acquisition,  guidance  with  video  detection  and  tracking  have  been  calculated  with 
the  projected  parameters.  For  tracking,  it  is  assumed  that  the  tracking  error  is 

CEP  0.35 
0  %  -  "  -  , 

m  \5W  Rim) 

and  a  pulse  integration  gain  of  17  dR  is  also  assumed. 

The  projected  performance  levels  for  acquisition  and  tracking  are  shown  in  Figures  lt> 
and  17  for  each  of  the  frequencies  and  a  1  kW  source.  Corresponding  curves  for  guidance 
were  not  plotted  because  power  adequate  for  the  other  two  functions  is  more  than  adequate 
for  guidance.  These  curves  show  some  improvement  over  those  of  Figure  15:  in  particular, 
they  show  that  it  is  possible  to  achieve  tracking  and  guidance  at  a  range  of  1  km  under 
more  stringent  conditions  than  before. 

The  calculations  in  the  study  show  that,  neglecting  angle  error  due  to  multipath  and 
glint,  94  GHz  is  the  best  frequency  for  a  multifunctional  tracking,  acquisition  and  beam- 
rider  guidance  system  using  state-of-the-art  performance  parameters.  For  a  system  domi¬ 
nated  by  multipath  effects,  either  140  or  220  GHz  appears  to  be  the  preferred  choice  fo: 
operating  frequency;  whereas,  without  the  presence  of  multipath  effects,  44  GHz  is  .«  better 
choice  for  development  of  first  generation  millimeter  guidance  systems.  This  result  is  to 
be  expected  for  performance  dominated  by  propagation  effects,  since  the  absorpt ion  due  to 
water  vapor  is  much  higher  at  140  and  220  GHz  than  at  94  GHz.  On  the  other  hand,  it  would 
seem  that  tracking  accuracy  would  favor  140  or  220  GHz  in  view  of  the  reduced  beam  dive: - 
gence  at  shorter  wavelengths.  For  a  conscan  tracker,  the  tracking  error  due  to  thermal 
noise  is  *  d^/vS/N,  and  the  \S/N  decreases  faster  than  t'^  as  the  wavelength  goes  to  2.14  ram 
(140  GHz)  and  1.36  mm  (220  GHz)  due  to  the  increase  in  mixer  noise  figure  and  atmospheric 
attenuation  coefficients  at  the  shorter  wavelengths.  Even  for  projections  based  on  extrapo¬ 
lation  of  available  source  power  and  improved  performance  of  quasi -optical  components,  the 
94  GHz  window  is  still  slightly  preferable  in  terms  of  tracking  accuracy  for  a  system  in 
which  multipath  is  not  a  consideration. 

On  the  other  hand,  if  the  multipath  model  adopted  is  valid,  there  are  indications  that 
94  GHz  will  not  provide  adequate  angle  tracking  error  over  typical  types  expected  to  be 
encountered  in  tactical  operations  in  adverse  environments.  If  the  trend  of  these  calcu¬ 
lations  is  consistent  with  projected  performance  from  realistic  terrain  models,  opeiation 
at  220  GHz  may  also  be  necessary  to  insure  an  adequate  accuracy  for  missile  impact  over 
the  desired  clear  weather  range  of  3  km.  If  this  is  the  case,  a  reduced  low  visibility 
operating  range  will  be  incurred.  Thus,  among  the  trade-offs  to  be  considered  will  bo 
that  of  adequate  clear  weather  operation  against  maximum  range  of  low  visibility  operation. 

A  94  GHz  system,  for  example,  might  perform  reasonably  well  over  a  1.5  -  2  km  range  interval 
in  a  wide  range  of  adverse  meteorological  and  terrain  conditions,  whereas  a  220  GHr  system 
might  extend  the  clear  weather  operating  range  to  2.5  -  l  km  while  restricting  the  adverse 
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weather  range  to  <  1.5  km. 


Figure  16.  Projected  performance  levels  for  Figure  17.  Projected  performance  levels 
acquisition  under  different  weather  for  tracking  under  different 

conditions  for  a  1  kw  source.  The  weather  conditions.  The  con- 

conditions  are:  moderate  rain  ditions  are  the  same  as  those 

4  mm/hr,  80°F;  heavy  fog,  30  m  given  in  Figure  16. 

visibility,  40'F. 

Preliminary  system  studies  have  shown  that  millimeter  guidance  systems  offer  a  poten¬ 
tially  significant  improvement  in  penetration  of  adverse  environments  which  limit  the 
visibility  of  electro-optical  guidance  systems.  Further,  there  is  the  potential  for  im¬ 
proved  accuracy  with  adequate  penetration  under  adverse  environmental  conditions  when  com¬ 
pared  to  guidance  systems  operating  in  the  microwave  frequency  band. 

The  state-of-the-art  will  currently  support  the  development  of  operational  breadboard 
beamrider  or  command  guidance  systems  at  94  or  140  GHz.  One  such  breadboard  system  has 
been  developed  at  BRL^O,  and  tests  performed  at  both  BRL  and  MIRADCOM  on  tracking  and 
guidance  link  operational  simulations  were  found  to  be  supportive  of  this  millimeter 
guidance  concept.  MIRADCOM  is  also  supporting  several  contractor-developed  millimeter 
guidance  systems  which  will  involve  development  and  testing  of  both  differential  guidance 
and  beamrider  concepts.  As  a  part  of  these  efforts,  the  evolutionary  development  of  cri¬ 
tical  subsystems  such  as  the  tracking  radar  and  guidance  link  will  provide  a  means  for 
performing  operational  tests  to  determine  multipath  effects,  measure  target  signatures, 
and  assess  various  schemes  which  have  been  proposed  for  missile  capture;  however,  a  very 
significant  concern  is  the  development  of  an  appropriate  target  acquisition  concept  tor 
an  autonomous  millimeter  guidance  system,  and  this  and  other  critical  technology  issues 
are  addressed  in  Table  9. 

A  basic  deficiency  in  the  current  millimeter  technology  base  is  the  lack  of  measured 
target,  terrain  and  atmospheric  data.  Without  these  data,  accurate  quantitative  estimates 
of  system  performance  in  realistic  tactical  environments  cannot  be  made  with  any  degree 
of  confidence.  Thus,  the  development  of  millimeter  measurement  systems  operating  in  the 
windows  at  94,  140  and  220  GHz  is  a  critical  stop  in  establishing  their  operational  limi¬ 
tations  in  terms  of  meteorological  conditions  and  reduced  visibility  situations  resulting 
from  smokes  and  other  battlefield-generated  aerosols.  The  acquisition  of  these  data  should 
be  a  major  part  of  a  program  to  develop  millimeter  guidance  systems. 
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ACT  IV I  Tt 
Multipath 


wavelength 

Selection 


Target 

Signature 


Target 

Acquisition 


Capt jre 


Acccleiat.tan 

Hardened 
Schottky  Diodes 

Near  m  wave 
Components 


Transmitters 


Target 

Tracking 

Techniques 

Clutter 

Characteristics 


Table  9.  Critical  Deamrider  Technology  Areas 


RATIONALE 

limits  tracking  accuracy;  Degrades  beam/ 
guidance 

Potential  of  Improved  tracking  vs  X;  source 
and  detector  performance  dependence,  multi¬ 
path  dependence;  optimum  search/acquisition 
wavelength 

Reduce  angle  error  and  signal  loss 


-ime  system  function 


Requ  rvd  for  Initiation  of  missile  guidance 


Required  for  missile  receiver 


Required  for  low  loss  sy'tem 


High  Power  Stable  Sources  Required 


Need  to  provide  high  accuracy  target 
coordinates  and  maintain  good  pointing 
accuracy 

Information  required  for  enhancement  of 
target-to-clutter  signals 


nw  x  radar/optical  tracker;  field  problems  - 
frequency  and  terrain  dependence 

Transmttcr/recei vers,  accurate  propaga¬ 
tion  measurements  in  smoke,  ram,  foo 
and  high  humidity;  supporting  meteorological 
measurements 

Transnltters/recei vers ,  Signal  Processing; 

Pul se-to-Pulse  target  reflectivity  data, 
function  of  tarqet  aspect 

Coherent  system  measurement s  to  feed 
computer  simulation  to  determine  signal 
processing  required,  e.g.  MT|,  vibrational 
signatures,  and 

Determine  required  beam  profile  vs  range 
for  capture,  near  field  measurement ;  deter¬ 
mine  concepts  for  multifunctional  capture/ 
track/guidance 

Trade-off  video/superheterodyne  for  missile 
receivers;  test  diodes  under  shock,  simulated 
flight  conditions 

Quasi -Optical  approach  -  circulators. 
Isolators,  modulators,  antenna  conscan 
schemes, diplexers,  rotary  joints,  mixers, 
local  oscillators 

EIO's  vs  Impatts;  develop  stability 
techniques,  chirp  techniques;  gyrotrons  - 
long  term  source 

High  precision  tracking  radar-spread  spectrum 
techniques  for  decorrelation  of  tarqet 
scattererv,  coherent,  frequency  agile  system 

Measurements  of  ffflW  clutter  characteristics 
terrain  masking  effects;  establish  dis¬ 
crimination  schemes 


Comments  on  other  systems  applications 


The  above  discussion  of  the  beamrider  indicates  some  of  the  factors  which  can  influ¬ 
ence  the  design  and  operation  of  a  NMMW  system.  The  same  type  of  considerations  must  be 
made  for  other  potential  NMMW  systems.  Only  brief  comments  can  be  made  here  on  some  of 
these  systems. 


Tarqet  designator  and  semi-active  seekers 


Target  designators  have  been  successfully  operated  at  1.06  pm  in  clear  atmosphere. 
However,  during  adverse  propagation  conditions  (mainly  fog  and  smoke),  the  IR  designators 
are  ineffective.  The  existing  IR  systems,  operating  at  1.06  um,  have  the  advantage  of 
having  all  laser  power  incident  upon  the  target.  This  is  important  since  the  reflectivity 
of  terrain  and  vegetation  exceeds  that  of  the  target  at  1.06  um.  For  NMM  wavelengths,  two 
techniques  have  been  investigated  [3] :  the  case  in  which  the  tarqet  fills  the  beam  and  the 
more  general  case  for  long  wavelengths  in  which  the  beam  size  exceeds  the  tarqet  dimen¬ 
sions.  Because  target  reflectivity  is  expected  to  exceed  background  reflectivity,  a 
largest  pulse  logic  is  postulated  for  designation  of  the  target  in  clutter. 


The  application  of  NMMW  command  guidance  to  anti-helicopter  or  anti-aircraft  missions 
appears  to  be  feasible.  Command  guidance  is  also  applicable  in  a  ground-to-ground  anti¬ 
armor  role,  but  no  clear  advantage  over  Beamrider  guidance  is  evident  at  this  time.  The 
technology  barriers  associated  with  this  guidance  concept  are  similar  to  those  for  the 
Beamrider,  and  the  general  conclusions  reached  for  the  beamrider  guidance  concept  apply 
to  command  guidance.  The  state-of-the-art  can  support  breadboard  system  development  at 
94  and  140  GHz,  and  additionally,  the  window  at  220  GHz  appears  to  be  applicable  if  source 
development  evolves  as  projected  by  current  technology. 


Target  surveillance  and  acquisition 


For  the  successful  operation  of  battlefield  NMMW  systems,  it  is  important  that  target 
acquisition  and  surveillance  techniques  be  applicable  to  NMMW  systems.  This  would  provide 
autonomous  operation.  For  the  optimum  weapon  system,  several  combinations  of  different 
wavelength  sub-systems  must  be  evaluated  to  obtain  the  optimum  combination  of  surveillance 
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acquisition,  and  guidance.  Among  the  systems  that  have  been  identified  as  potential  NMMW 
applications  are: 

battlefield  ground-to-ground  target  acquisition; 

airborne  surveillance  by  RPV  radars  of  land-based  targets; 

long-range  surveillance  and  target  acquisition; 

target  acquisition  for  NMMW  airborne  missile  seeker  systems; 

horizon  search  radars; 

shipboard  surveillance  radars; 

track-while-scan  systems  for  short-range  self-defense  systems 
(point  defense) ; 

restricted  scan  target  acquisition  for  missile  systems. 

A  large  number  of  equipment  and  technique  developments  are  necessary  to  ascertain  the 
applicability  of  tPlM  waves.  Techniques  for  search,  recognition  and  classification  have 
been  identified  [3]  but  more  detailed  investigations  are  needed  to  implement  these  tech¬ 
niques.  Greater  transmitter  power,  improved  receiver  sensitivity,  new  low-loss  components 
and  light-weight  rapid-scan  antennas  are  necessary.  Source  stability  for  coherent  oper¬ 
ation  must  be  greatly  improved.  For  implementation  of  the  recognition/classification 
techniques,  target  characteristics,  clutter  effects  and  atmospheric  effects  are  among  the 
phenomenology  that  must  be  thoroughly  documented. 

Boost  phase  plume  detection 

Calculations  have  been  performed  to  determine  the  feasibility  of  detecting  missile 
plumes  during  the  boost  phase  of  the  vehicle  [22].  In  the  altitude  regime  of  30-100  km, 
it  has  been  shown  that,  for  a  solid  propellant  system,  molecular  species,  e.  g.  H20  and 
H  Cl,  emit  to  permit  radiometric  detection  from  aircraft  or  satellites.  Airborne  radio- 
metric  observations  are  required  to  confirm  predictions.  The  necessary  passive  technology 
is  developing  but  must  be  extended  to  wavelengths  as  short  as  -350  um.  Independent 
analysis  of  the  bus  stage  indicate  that  spectral  line  opacity  is  sufficiently  high  to  be 
detectable  in  occultation  against  the  earth's  radiation. 

Conclusions 

From  the  studies  which  have  been  performed,  it  can  be  concluded  that  the  NMMW  region 
offers  a  compromise  for  good  resolution  under  adverse  propagation  conditions.  Most  sys¬ 
tems,  which  have  been  investigated,  profit  from  the  narrow  antenna  beams  available  at  NMM 
wavelengths,  but  are  limited  by  atmospheric  conditions.  In  some  cases,  the  range  limita¬ 
tions  imposed  by  the  atmosphere  serve  as  an  advantage  for  convert  operation. 

In  order  to  utilise  the  NMMW  spectral  reeion  properly,  considerable  technology  must 
be  developed.  Hioh  -power  sources,  low  loss  components,  precision  antennas,  sensitive 
receivers  and  stable  local  oscillators  are  priority  devices  to  allow  systems  operations 
to  be  performed  at  TV'  wavelengths .  For  man'  systems,  highly  coherent  transmitters  are 
required.  This  necessitates  development  of  phase-locking  and  injection-locking  technology 
for  high-power  sources.  A  driving  force  in  extending  military  operation  to  NMM  wavelengths 
is  the  prospect  of  sr.all-size  systems  where  space  and  antenna  apertures  are  limited.  Low 
cost,  which  is  not  a  current  characteristic  of  NMMW  components,  is  expected  to  be  an  ulti¬ 
mate  achievable  goal. 

One  cannot  exoect  to  achieve  everything  that  a  system  demands  by  employing  NMM  waves, 
and  those  working  in  the  field  are  not  doinc  this.  Phenomenology  in  the  form  of  atmos¬ 
pheric,  terrain/clutte - ,  target,  and  materials  characteristics  must  be  thoroughly  developed 
for  comparison  with  operation  at  other  wave lenc ths .  Thus  far,  studies  have  estimated 
certain  characteristics  and  projected  reasonable  operational  parameters  for  devices.  From 
the  considerations  that  have  been  made,  optinur.  operation  can  be  expected  from  hybrid 
IR/NMMW  systems  which  will  utilize  the  best  of  both  spectral  regions. 

This  work  has  been  sponsored  in  part  b’  the  Harry  Diamond  Laboratories  through  the 
Army  Research  Office  Contract  No.  DAAG29-77  -C-1326. 
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